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Introduction 


___  Chapter  1. 

i  beery  ^  tn^rqy 

'  Introduction 


i  f(  vZ/  *“JU/  C  t  »*y>du 


(r 


*>  This  investigation  has  the  purpose  of  performing 


seif-consistent  energy  band  calculations  on  some  of  the  II-VI 
compounds,  such  as  cadmium  sulfide,  zinc  oxide,  and  sine  sulfide. 


There  are  many  practical  and  theoretical  reasons 


fer  interest  in  these  materials. 


'-jThese  materials  have  applications  as  phosphors  (ZnS 

and  CdS),  as  infrared  detectors  (ZnS),  in  photovoltaic  cells  (CdS- 
n 

Cu  S  or  CdS-CdTe  hetero junctions) ,  in  batteries  (ZnO),  as  FETs 
2 

(CdS),  In  heterojunction  lasers  (CdS),  and  even  as  acoustic 
amolifiers  (CdS).  ^  ,  1>C$  )  *— 

J 

TJhiey  arc  of  considerable  theoretical  interest  as 


well,  as  examples  of  wide  band  gap  solids,  intermedia tc  between 


the  covalent  III-V  semiconductors  and  the  more  ionic  I-VIJ 


insulators.  They  serve  as  a  test  of  calculational  methods  previously 


used  on  rhe  lighter  III-V  compounds  and  other  compounds  simpler 


than  these  II-VI  compounds. 


3 


In  particular,  the  occupied  el-levels,  which  lie  in  energy  near  the 
top  of  the  valence  levels,  complicate  the  picture.  The  positions  of 
the  d-bands  have  not  been  correctly  predicted  in  previous  ab  initio 
calculations,  and  in  fact  some  methods  fall  qualitatively  by  not 
obtaining  the  proper  ordering  of  the  valence  s,p  and  d  levels. 

The  II-VI  compounds  have  also  been  of  theoretical  Interest 
os  excitonic  systems.  High  densities  of  such  electron-hole  complexes 
have  been  produced  and  observed  in  CdS.  The  electron-hole  liquid  has 
been  observed  and  studied  In  CdS  as  well. 

In  addition,  some  very  interesting  magnetic  properties 
in  CuCl  and  CdS  have  recently  been  observed.  This  is  the  background 
of  these  experimental  studies.  Anomalous  diamagnetism  has  been 
observed  in  CuCl  under  rapid  pressure  and  temperature  cycling.  In 
some  samples,  up  to  50£>  of  the  magnetic  flux  is  excluded.  This  occurs 

to 

at  high  temperatures  (150  k!  )  and  cannot  be  explained  *y  ordinary 
diamagnetism.  The  only  currently  known  phenomenon  that  can  explain 


this  is  superconductivity. 


9 


It ‘a  clear  that  this  would  have  to  be  superconductivity 

of  different  origin.  Acoustic  phonons,  the  known  mechanism, 

are  too  low  in  energy  to  bind  electrons  at  such  high  temperatures. 

Models  using  some  other  form  of  interaction  to  pair  electrons 

seem  to  be  necessary.  Some  form  of  excl tonic  binding  is  proposed 

in  many  of  these  models.  One  of  the  most  well-known  such  models, 

the  Allender-Bray-Eardeen  model,  requires  a  metal  and  a 

semiconductor  to  exist  in  very  close  proximity.  Since  CuCl  can 

disproportionate  via  the  reaction  2CuCl  Cu  +  CuCl  ,  this  model 

2 

has  appealed  to  many.  Kovr/cr,  cadmium  sulfide,  which  has  been 
observed  to  have  only  one  valence  state  and  is  net  believed  to 
disproportionate,  hs.s  shown  similar  anomalous  diamagnetism.  The 
law  of  parsimony  is,  therefore,  against  the  metal-scmiconductor 
sandwich  idea.  Also,  in  both  CuCl  and  CdS,  impurities  clearly 
play  a  major  role.  Sufficiently  clean  samples  display  no 
interesting  behaviour  at  all.  And  to  make  life  even  more 
interesting,  cadmium  sulfide  also  becomes  a  ferromagnet  under  a 


high  (40  kilogauss)  applied  field. 


Another  model,  proposed  by  Bishop  end  Overhauser  (  ),  Involves 

the  interaction  mediated  by  the  optical  phonons.  Thin  mechanism 
has  deep  attractive  potential  wells,  about  10  mev  in  depth  for 
a  spacing  of  250  angstroms.  This  suggests  the  possibility  of 
bound  states,  and  the  pairing  is  considerably  stronger,  around 
a  factor  of  ten,  than  the  pairing  from  acoustic  phonons  in 
the  typical  Cooper  pair  in  the  highest  temperature  superconductors 
known ,  at  approximately  20  degrees  Kelvin.  This  idea  is  even  more 
intriguing  in  the  light  of  the  fact  that  both  CdS  and  CuCl  are 
strongly  polar  compounds.  This  suggest  that  these  phenomena 
would  r.ot  be  seen  in  the  less  polar  IIX-V  compounds  and  in 
elemental  semiconductors  such  as  Ge  and  Si.  They've  been 
investigated  most  thoroughly  and  indeed,  such  effects  have  not 


been  observed 


Chapter  One 


The  calculations  discussed  here  are  based  on 

Hartree-Focic  theory,  the  most  common  approximation  to 

the  exact  non-relativistic  theory  of  a  many-elcctron  system. 

I'll  discuss  the  theory  and  its  applications  to  crystalline 

solids.  The  specialization  of  Kartree-Fock  theory  used 

in  these  calculations,  the  method  of  local  orbitals, 

deserves  ana  gets  a  chapter  of  its  own. 

The  original  problem  is  finding  exact  solutions 

to  the  Schrodinger  equtation  for  a  many-electron  system. 

This  is  impossible,  if  the  universe  is  as  small  as  we  think 

it  is.  Consider  the  wavef unction,  which  is  a  3N  dimensional 

function  if  there  are  N  electrons.  Dividing  each  axis  into, 

say,  100  units,  in  order  to  numerically  integrate  and 

6n 

differentiate,  we  need  to  record  10  entries.  The  same 


problem  arises  if  we  try  to  describe  the  wavefunction  by 


sets  of  orthogonal  functions.  10  is  just  too  big  a  number. 

Quantum  mechanics  is  agaiar.st  us.  The  upper  limit  for  the 

2 

rate  of  information  flow  in  a  computer  of  mass  K  is  Me  /h 
(h  ss  plnnck's  constant,  c  =  speed  of  light),  in  bits  per 
second.  It's  not  hard  to  see  that  macroscopic  physical 
systems  cannot  b.e  exactly  simulated  by  digital  computers. 

I  emphasise  the  v/ord  digital,  since  the  electron  is  an 
excellent  analog  of  Itself. 

Such  a  digital  simulation  would  be  useless  even 
if  it  could  be  done.  As  suggested,  the  real  solid  simulates 
Itself  perfectly.  The  goal  of  physics  is  to  explain  some 
set  of  occurrences  in  terms  of  a  few,  relatively  simple 
ideas,  rather  than  making  a  model  isomorphic  with  the  universe 
end  just  as  confusing. 

* 

We  begin  to  clear  away  some  complications  by 

using  the  Eorn-Oppenheimer  approximation.  This  neglects 

any  relationship  betvreen  the  motion  of  the  electrons  and 

the  motion  of  the  nuclei.  This  a  valid  simplification,  since 

3  5 

the  nuclei  range  from  10  to  10  times  heavier  than  an  electron. 


We  can  express  the  K-electron  wave function 


as  a  function  of  the  electronic  coordinates  alone? 


f  =  tc*,*  ...  *  ) 

K  f  * 


(t.l) 


Here  contains  the  position  and  spin  coordinates  of  the 


1th  electron.  We  write  the  Schrodinger  equation  for  this 


Kftvef unction: 

^  =  ^  % (*<* (l,2) 

Here,  H  is  the  non-rclatlvistlc  Hamiltonian  expressed  in 


atomic  units.  In  this  system  of  units,  a  1,  e  ~  2,  and 


the  mass  of  the  electron  is  one-half.  The  unit  of  distance 


is  one  bohr  (.52 9  angstroms)  and  the  unit  of  energy  is  the 


rydberg  (13*6  cv). 


H  ’  H,  ’K 

K*  “I  Jr  f  -l  2  2  rJr-rg, 

2  1  J  J  1  J  lH-fU 

U  jV  I  il 

U  -  .  r  1>J  U 

H;+~  Z  -P- 

1 2  x  Of-  S  ) 

J  * 


(1.3) 


(1.4) 


H  is  independent  of  the  electron  coordinates,  and  has  been 
n 

assumed  to  be  a  constant,  as  per  the  Eorn-Oprenheimer 


approximation. 


c; 


The  upper  case  characters  refer  to  nuclear 
properties:  2  is  the  atonic  number,  R  is  the  nuclear  position 
(  usually  assumed  to  ce  fixed,  like  Pimlico)  and  M  is  the 
nuclear  mass.  The  lower  ease  characters  describe  properties 
of  the  electrons.  y( , is  the  coordinate  of  the  ith  electron. 

A 

The  atonic  numbers  of  the  atoms  under  consideration 
are  relatively  low,  and  since  in  any  event  it  is  the  valence 
electrons  that  arc  of  primary  interest,  the  non-relativistic 
nature  cf  the  Hamiltonian  may  be  acceptable. 

The  fundamental  approximation,  nov;  applied,  is  the 
independent  particle  model.  In  truth,  in  the  real  solution, 
the  variables  are  not  separable,  l/e  must  assume  that  they  are, 
or  more  exactly  that  the  true  wavefunctlon  can  be  well- 
approximated  by  this  model.  The  independent  particle  model  “says 
that  the  electron  is  acted  on  by  the  average  of  the  other 
electrons.  The  N-electron  wavefunctlon  is  expressed  as  a  product 
of  one-electron  wavefunctions,  or  as  a  linear  combination  of 
such  products.  A  simple  example  is  the  Hartree  wavefunction: 


/O 


Since  electrons  are  fermions,  this  K-electron  wovefunctlon  must 
be  nntisymroetrlc.  Any  interchange  of  two  particles  must  reverse 
the  sign  of  the  wavefunction.  The  simplest  way  to  antisymmetrize 
is  to  use  the  determinant  of  a  matrix  whose  elements  are 

( «  •/ 

' '  ' 

Not  all  such  are  spin  eigenstates,  but  fortunately  we  liave 
A* 

here  closed-shell  systems  that  can  be  expressed  as  single- 


detcrwlner.tal  wavef unctions. 

For  any  arbitrary  ^  ,  we  can  find  the  expectation 
value  of  the  Hamiltonian.  It  is  a  sum  of  one-  and  tv/o-  particle 


Integrals : 


P  is  a  permutation  operator,  constructed  so  that 

o  * 

/  I  a  Zy  “  /  >  A  y  .  The  terms  i=j  are  not  a  oroblem  since  the 

/,/  V  J  ■> 


self-energy  cancels  the  ’self-exchange*. 


The  expectation  value  of  the  Kami  1  ton ian, <  /  P  f  /,  /,  can  be 

•  »(  /V 

i 

minimized  by  varying  the  spinorbltals.  It  has  been  shown  that 

< l  Hi  t>z  E  a -9) 

<  t  i  t) ' 

where  E  is  the  energy  of  the  exact  ground  state. 

We  vary  the  orbitals,  but  take  care,  to  keep  them 
orthogonal,  for  convenience  and  without  loss  of  generality.  This 
is  done  by  using  Lagrange  multipliers.  The  objective  is  to  obtain 
the  lowest  upper  bound  to  the  exact  non-relativistic  energy. 


We  define  the  functional 


and  require  that  the  variation  of  L  =  0  for  1  =  1,2,...N  . 


This  yields  a  set  of  K  coupled  nonlinear 


integrodifferentinl  equations  known  as  the  Kartrce-Fock  equations. 

{•' -1  /Ipr  +  i  CJ}  -  *; )}  % <x)  =fv^(1-u> 

J  and  K  have  the  familiar  Coulomb  and  exchange  operator  forms. 

J  J  2 

J.Aw  =  f  f  Ip^lL  S/}  0.to  (1.12) 

K;  &-M  ctfb')0£V)  yyj  O.M  (1.13)  ■ 

Iprr  "  ) 

These  equations  may  be  written  in  matrix  form; 


(1.1/0 


n 


F  is  the  I!xN  matrix  of  the  Fock  operator.  For  closed 
shells,  it  is  h'ermitian,  so  S  ray  be  diagonalized  by  a  unitary 
transformation.  Fortunately,  the  II-VI  compounds  have  closed 
shells.  We  aay  write,  F //>  =*  £ ,  / 1  >  .  (1.15) 

A 

Combining  this  expression  with  equation  1.7  gives  the  Hartree-Fock 
ground  state  energy, 


F. 


N 

KF 


(i/2)l  f  £;  +  </(//.<>}  .  (i.i6) 


The  are  the  single-particle  energies,  almost.  Koopman's 


A 

theorem  states  that  given  a  variationally  stationary  state 

formed  from  H  simultaneous  spinorbital  (pseudo) eigenfunctions 

of  the  K  particle  Fock  operator  F,  the  states  9  and  9  , 

A/*-/  fu-j 

formed  by  altering  the  list  of  occupied  eigenfunctions  of  F  by 
one  entry,  are  stationary  with  respect  to  further  variation  of 
that  spinorbital  whose  occupancy  has  been  altered.  If  we 
neglect  relaxation  effects,  £*  is  the  entire  change  in  system 
energy  when  an  electron  is  added  or  removed. 

In  order  to  equate  the  C .  5  with  single-particle 

A. 

energies,  we  must  also  assume  that  relaxation  effects  are  small  - 
that  the  other  occupied  orbitals  change  only  slightly  when  an 


n 


There  are  such  changes,  since  the  Pock  operator  is  a  functional 
of  nil  of  Its  own  solutions,  and  changing  one  of  then  must 
change  all*  The  question  is,  how  much?  These  relaxation  effects 
will  be  discussed  in  a  later  section. 


The  Method  of  Local  Orbitals 

The  method  of  local  orbitals  Is  a  variational 

technique  In  which  solutions  to  the  Kartree-Fock  equations 
* 

are  sought  for  a  small  subsystem  of  the  system  of  interest, 

with  the  caveat  that  all  such  solutions  found  must  be  contained 

within  the  occupied  Fock  space  of  the  original  system.  Sets  of 

such  local  orbitals,  which  arc  not  orthogonal,  can  then 

be  assigned  to  each  subunit  of  the  larger  system  using  the 

translational  invariance  of  the  Pock-Dirac  density  matrix, 

thus  spanning  the  occupied  Hilbert  space  of  that  system. 

A  single  noniterative  rotation  of  the  local  orbitals  within 

this  occupied  space  then  yields  the  exact  self-consistent 

solutions  for  the  original  Pock  operator. 

If  £  0A-  .  A  z  Ijl’j"'  h  is  the  set  of  local 

orbitals  satisfying  the  Kartree-Fock  equations  for  an  n-clectron 

system,  then  F  0:  ~  jl  /I.  ,  &•  ■  ,0.1) 

A  J  ‘ J 


The  Fock-Dirac  density  o aerator  for  the  equation  is 


/  J 


(3.3) 


•«/  f' 

-  T  &■  S-.  b.  (3.2) 

-1  V  * 

S  is  the  Inverse  of  the  overlap  matrix  s. 

S,  -  j  &jx)  (?'(x)  c/ k  J3.3) 

*  y  A  / 

2  ^ 

Rho  (j*)  is  a  projection  operator  (j>  *j>  ) 

chosen  such  thatj?&*  Q  for  any  orbital  ^'in  the  occupied  space 
andjv/ifaO  for  any  orbital  %  in  the  virtual  space.  This  means 
that  for  any  one-electron  operator  L,  the  projection  of  L 


onto  the  Knrtree-Fock  manifold  will  satisfy 


y  L> 6;  *£,  y-j  % 

for  any  occupied  orbital^',  and^,  Lj,  f=0  for  any  virtual  orbital  5* 


Another  approach,  rigorously  developed  by  Gilbert  (  ), 

begins  by  introducing  modified  iiartree-Fock  equations  of  the  form 

(  f  +jhj>  &'  (3.5) 

Above,  A  is  an  as  yet  unspecified  Hermltian  operator  and  rho 
is  the  ’diagonal '  (  x  equals  x  prime)  form  of  the  Fock-Dirac 
density  operator  (  ).  We  see  that  rho  is  idempotent,  Hermltian, 

and  projects  onto  occupied  Pock  space.  Because  of  this  last 
property,  it  can  be  shown  that  the  occupied  eigenfunctions  of 


the  modified  Har tree-rock  equation  lie  entirely  within  the 
occupied  Fock  space  of  the  unmodified  Hartree-Foek  equations, 

,  /  C  $j  f  A+- 

SO,  Ji  (  Q  ft Jff  (3»^) 

There  is  no  Koopman's  theorem  for  this  modified 

« 

equation,  and  so  the  f.s  cannot  be  interpreted  as  being  nearly 
the  particle  energies,  as  they  can  in  the  regular  Hartree-Fock 
equations.  The  total  energy  of  the  subsystem  still  has  significance 
since  a  unitary  transformation  within  occupied  Fock  space 
leaves  the  trace  of  Hermitian  operators  such  os  the  Hamiltonian 
invar lent. 


The  properties  of  the  ^*will  vary  with  the  choice 

a 

of  A.  We  are  looking  for  functions  that  are  localized  on  specific 
lattice  sites  of  a  crystal.  In  order  to  decide  which  operator 
best  fulfills  such  a  purpose,  we  need  to  consider  the 
minimization  of  a  functional  . ... 

subject  to  the  variational  constraints  <  /lj>  -i.il  =0; 


P 

J 


U)  ~JJ)  -I'j, 


=0 


(3.7) 


Here  we  can  use  the  method  of  Lagrange  multipliers.  The  Integral 


constraints  of  equation  3.7  arc  taken  into  account  by  the 


rarameters  >?■*  *  and  the  constraints  that  involve  explicit 

'r 

functions  of  X  In  equation  3.7  are  taken  care  of  by  I-agronge 


r.  o  < 

functions  /£  (x)  and  ,^(x). 


We  obtain 


?{  -  f  lj<  **  ip  -  i .  i  } 

*  1  Vy‘*l  -/•  '  V 

■"  ,2.  '^{x)  Uxlj?  —  <Z  /  )  ’  ~Z  ^  / *2  (j>ts  >  ~  J  ^.1 

=  /’  M-lg  1  ~  J  2  .!*>  ^  ^  o;?  *  ^  -o 

x.  ;=/  /  /<  J  ** 

Kultiplying  on  the  left  by  (1 —p)  •  we  see 


(l-jrt  l^>  &0 


and  T  J>  -*0  for  all  ^ 

and  so  +  ^(x)  -J>f2  J&to}  ~'4j:c) 

Substituting  equation  3.H  into  3.3.  we  get 

y/j&ta]  “  I  7  //> 

J  *  J-! 

which  may  be  rewritten  os  Q)ky  =  y  2.  I  -y } 

&  MJ  ' 

where  £>  =  r>  cf.J))  <'  <>/  = 

JVa#  '  v 

/:  unitary  transformation  gives 


(3.9) 


(3.10) 


(3.H) 


(3.12) 


(3.13) 


(3.14) 


6  //<•>  =  £'/«'> 
A 


(3.15) 


Since  we  now  know  how  to  minimize  any  chosen  Kennitlan 


functional  *<0-  ,  we  can  now  get  lack  to  tha  true  problees.  Let 


A  -  G  •  F.  *o  that  [  (  F  p  G2  ~  O  (3.16) 

We  cen  now  rainlmlzo^jf^Jand  solve  equation  3«5  at  the  same  time, 
and  so  may  picU*H<^to  satisfy  physical  reasoning. 

We  choose  A  so  that  the  solutions  are  localized  on  a  site, 


usually  by  picking  A  to  be  some  potential  well  centered  on  the 


site.  The  rock  operator  can  be  separated  into  P  ,  the  atomic 

s 

Pock  operator  for  the  site  s,  and  U  ,  an  external  potential 

s 

operator  for  the  site  s.  We  usually  pick  the  localizing  potential, 


A,  to  be  -U  .  The  local  orbitals  equation  is  nov; 
s 


£  Fi  +  ^  -  f.  (3.17) 

This  isn't  an  obvious  improvement  over  the  original 
Hartree-Fock  equation,  but  it  allows  a  systematic  approximation 
that  simplifies  greatly  the  Hartree-Fock  problem  without  inducing 

to 

undue  error.  The  left  side  of  equation  347  can  be  analyzed  by 
the  order  of  the  intersite  overlap.  The  Fock-Dirac  density 

y. 

operator  is  p  =  [  [  9  V  ^  Q  (3.18) 

J  Si  7}  5'  ^ 

with  V  the  Interatomic  overlap  matrix. 


Ja 


s 


*4 


As  discussed  in  Lowdln  {  ) ,  \J  _iAay  be  expanded  in  a  power 


-/ 


series  in  V  _  ?  V  ■  ]  j  -  V 


Vv 


(3.20) 


,  /  ><*«'■! 
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Expanding  the  local  orbitals  equation  and  discarding  second 


order  and  greater  in  intersitc  overlap  ( 


)  gives 


U 7 


(3.21) 


As  self-consistency  is  approached,  the  last  two  terras  of  equation 
3.21  tend  to  cancel.  Therefore,  the  Pock  space  rotation  changes 
Bloch  v/avefunctions  into  atcaic-like  wavefunctlons  that  are 
eigenfunctions  of  the  atomic  Pock  operator. 

We  procede  to  localize  a  physically  realistic 
number  of  electrons  on  each  site.  In  cadmium  sulfide,  we  choose 

m 

to  have  46  electrons  on  the  cadmium  site  and  38  on  the  sulfur, 
approximating  an  ionic  charge  distribution  (  overlap  and  diffuse 
valence  orbitals  permit  complete  covalency  if  needed,  however). 

V,'e  have  achieved  several  advantages.  The  problem  is 


computationally  simplified,,  since  the  equation  being  iterated 


to 

does  not  have  to  portray  the  full  complexities  of  the  solid, 
and  since  the  local  orbitals  are  recognizably  modified  atoalc 
orbitals  and  aid  in  picturing  chemical  bonding.  This  approach 
has  use  as  a  beginning  in  looking  at  amorphous  substances,  which 
have  only  a  short-range  order. 

The  local  orbitals  are  calculated  using  the  program 
called  LOPAS  written  by  Kunz.  A  basis  function  expansion  of 
Slater  orbitals  (STOs)  is  used  for  the  radial  part  of  the  local 


orbital  while  a  spherical  harmonic  describes  the  angular  behavior. 

Bt  (?)  -  TO  &  (3.22) 

We  usually  take  the  /u  and  2^  .  from  Eagus  et  al  (  )  where 

V 

optimized  basis  sets  have  been  computed  for  the  atomic  or  ionic 
system  appropriate  to  the  solid  in  question.  In  many  cases,  we 


will  change  this  basis  set  to  achieve  localization,  considering 


at  the  same  time  the  total  energy  of  the  subsystem. 

The  expansion  coefficient  is  solved  by  the  matrix 

/ 

method  of  Roothean  (  ). 


Equation  3*21  is  solved  repeatedly  until  the  self- 


consistency  is  below  some  predetermined  value,  normally  one  part 

-6 

in  10  for  the  charge  density. 


a’2  have  only  weed  the  long-range  pars  of  the  external  potential 


as  the  localising  potential.  So,  for  ionic  substances,  the 
electrostatic  potential  will  be  screened  while  the  short-range 
effects  of  the  inner  shells  will  change  the  local  orbitals,  Some 
anions  will  not  localise  in  this  situation  since  it  is  the 
Kadclung  part  of  ‘the  potential  that  stabilises  them.  Since  A  is 
essentially  arbitrary,  v<t  are  free  to  niter  it  in  whatever  way 
that  achieves  optimal  localization,  and  we  often  use— ^  jy 
in  equation  3.21  instead  of  tjx 


Energy  Band  Theory 


The  Hertree-Fock  energy  hands  ore  obtained  by  plotting 
the  one-electron  energies  derived  In  equation  1.16  against  the 
crystal  momentum  k.  Our  calculations  make  use  of  the  symmetries 
of  the  crystal.  In  the  case  of  these  calculations,  it  is  the 
sine-blende  crystal  structure. 

Bloch's  theorem  requires  that  the  vmvefunctlons  of  the 
solid  satisfy 


,u  _  ^  V 

I  “  -a-  AX  C  v )  (2.1 ) 

hk 

xl  (  y  r  R  )  - 

^  k 
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hk 
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MS 


where 

for  any  lattice  translation  vector  R.  The  crystal  momentum  is 


again  k  and  n  is  the  band  index.  V/avefunctions  that  satisfy 


this  condition  are  Bloch  wavef unctions. 


If  v;e  look  at  the  set  of  R  ,  we  can  see  that  there 


n 

is  an  element  of  R  that  1ms  a  minimum  length,  corresponding  to 

n 

the  fundamental  lattice  spacing.  Considering  this,  and  the  fact 

that  R  is  o  group  under  translation,  we  can  sec  that  there  is 
n 

a  max inum  k  when  we  transform  into  momentum  space. 


Th*5rofor*.,  fcher*  i«  ».  finite  vision  of  non-'vju  lvale.ee  v-noJnts  - 


the  first  Erlllouin  xone.This  may  be  reduced  to  on  '  irreducible 
wedge'  by  symmetry  considerations. 

The  solution  to  the  Knrtrce-Fock  equations  is  lengthy 
for  crystals,  and  the  successive  approximation  nature  of  the 
calculations  requires  that  the  equation  be  solved  several  times. 


The  Kar trce-rock  equation  can  bo  rewritten  in  terms  of  the 


Fock-Dlrac  density  ratrix(  )t 

I  *  n£r,  f  ><*']  &  ») 


where  p(Jjjt')-  I  )&*&') 


(2.2) 


(2.3) 


If  a  self-consistent  density  matrix  is  available,  then 
equation  2,2  r.eed  only  be  solved  once  to  obtain  the  eigenvalues 
and  eigenfunctions,  Adams  (  )  and  Gilbcrt(  )  have  shown  that 

the  density  matrix  (x,x*)  Is  the  same  for  equations  2.2  and  >.21  t 


-/  ^  t 


v  t/.r)  *  Xsj )  e*w)  *  1 1  £  V '  6 

/  3  <  •  t,  j  j*r.  r„ 

*  )  3  7  * 


(2.4) 


So  ;;c  have  the  required  density  matrix,  and  can  obtain 


crystal  wavefv.nctlons  from  the  already  obtained  local  orbitals 


by  constructing  the  Fock  operator  and  solving  equation  3.2i  Just 


once.  The  matrix  elements  of  F  ore  calculated  to  first  order  in 


interatomic  overlap  ,  consistent  with  the  calculation  of  the 


local  orbitals. 

The  Bloch  functions  $_(r)  arc  expanded  in  a  basis  set 

KJ< 

in  which  the  basis  functions  have  the  form 


This  linear  combination  of  basis  functions  technique  is  similar 
to  the  famous  linear  combination  of  atomic  orbitals  technique 
(LCAO)  except  for  the  fact  that  the  are  not  free-atom  orbitals, 
but  are  the  local  orbitals  obtained  previously.  Note,  the  local 
orbitals  are  occasionolly  enriched  for  the  case  of  virtual  states. 
The  LOFAS  program  will  not  indicate  basis  functions  appropriate 
for  a  virtual  state  if  that  virtual  state  is  of  a  different 
angular  momentum  type  than  the  occupied  states?  for  instance, 
a  case  in  which  the  occupied  states  are  all  s  and  p-like,  while 
the  first  virtual  state  is  d-llke.  In  this  case,  basis  states  * 
for  the  virtuals  are  added,  consisting  of  spherical  harmonics 
multiplied  by  single  STO*s  that  are  chosen  tc  have  small  overlap. 

Since  the  crystal  momentum  k  is  a  good  quantum  number, 
the  Fock  space  is  diagonalized  into  separate  spaces  for  each 
reciprocal  lattice  vector.  The  integrals  necessary  to  perform 


I > 

the  calculation  are  l:-\ndep^.  :*  «*t  and  need  only  r er formed 
once.  Hulticenter  integration  is  done  by  the  Lows  in  o(-fur.ction 
expansion  method  (  ),  One  of  the  sites  is  chosen  as  center  and 

all  functions  are  expanded  in  terms  of  spherical  harmonics  centered 


on  that  site.  The  calculations  are  performed  using  the  programs 


K3AKD  and  K30KE  uritten  by  Kunz.  The  output  consists  of  energies 


and  coefficients  of  the  basis  functions  at  20  selected  points 


of  the  Brlllouin  zone. 


•  Correlation  Corrections 


As  previously  mentioned,  Hartree-Fock  theory  does 
not  give  exact  solutions  to  the  true  many-electron  Kamil tonlan. 
Since  we  are  using  a  single  determinental  wavefunction,  electrons 
are  affected  only  by  the  mean  field.  Electrons  of  the  came  spin 
have  some  of  their  true  pair  interactions  taken  Into  account, 
but  there  is  no  pair  correlation  at  all  btween  two  electrons 
of  opposite  spin.  It  is  obvious  that  the  Coulomb  force  between 
two  electrons  is  independent  of  spin,  and  should  keep  electrons 
from  the  near  neighborhood  of  any  other  electron. 

We  say  that  the  motion  of  electrons  of  opposite 
spins  is  uncorrelated,  and  that  electrons  with  like  spins  are 
Incompletely  correlated.  The  correlation  energy  is  usually 
defined  as  the  difference  between  the  energy  obtained  from 
our  mean-field  Hartree-Fock  theory  and  the  exact  non-relativlstlc 
energy  of  the  system.  This  exact  energy  could  in  principle  be 
determined  by  using  a  computational  technique  called  configuration 


interaction 


In  this  method,  a  wave function  is  used  that  is  a  linear 
combination  of  determinants.  It's  impractical  for  all  but  the 
smallest  molecular  and  atomic  systems,  and  is  for  too  time- 
consuming  to  be  useful  in  solids. 

For  solids  in  general,  ond  specifically  in  the  case 
of  the  II-VI  compounds  under  discussion,  we  attempt  to  estimate 
the  correlation  correction  to  the  Kartrce-Fock  calculation.  The 
correlation  correction  for  the  cnc-electron  eigenvalue  is  the 
error  from  Koopmaru  '  theorem  (  )  end  from  the  independent 

particle  model. 

According  to  Koopmans*  theorem,  an  orbital  eigenvalue 
is  approximately  the  difference  beyween  the  Hartree-Fock  energy 
of  the  N-electron  system  with  the  level  occupied  and  the  Kartree- 
Fcck  energy  of  the  K-l  electron  system  with  that  same  level 


unoccupied.  These  Kartree-Fock  excitation  energies  are 


0  r-W  ~  <"-0 

l 


(4.1) 


for  states  that  are  occupied  in  the  Kartree-Fock  ground  state 


f  Afrt)  (A/) 

f  -  t  -  F 
\,-u'  Jt  *  //  /- 


(4.2! 


for  virtual  states  of  the  Hartree-Fock  ground  state. 


When  we  attempt  to  improve  the  Hartree-Fock 


Approximation,  we  wish  to  keep  the  goneral  hand  scheme,  so 

(4.3) 

=  f  <*>”>.  F«)  (4.4) 

c  •  "" 


_  (*)  _  t  At—) 
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where  the  eigenvalues  are  the  exact  eigenvalues  of  the  system. 

We  assume  that  the  correlation  energies  are  o  small  perturbation 

(L) 


of  the  original  Hartree-Fock  bands,  and  we  can  then  write  E 
<L) 


as 


(L)  (L) 

E  +  E 


(4.5) 


HF  c 

(L) 

here,  E  is  the  total  correlation  energy  of  a  system  with  L 


electrons.  Equations  4.3  end  4.4  may  now  be  written  as 


0  (K) 

OM.) 

(4.6) 

e  =  e  4*  (E  - 

E  ) 

nk  nk  c 
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0  (N) 

(K-l ) 
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e  a  e  +  (E 

E  ) 

nk  nk  c  c 


A  useful  approximation,  valid  for  nonmetals  with 

valence  band  width  lrss  than  the  Hartree-Fock  optical  band  gap, 

«• 

was  developed  by  Fantelides  et  al  (  ).  They  showed  that  equations 

4.6  and  4.7  may  be  replaced  by 


0  (N-J  ) 

E  =  E  +  E  (h) 

nk  nk  nk 

(4.8) 

0  (K) 

E  a  E  +  E  (e) 

rik  nk'  rik' 

(4.9) 
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CO 

In  these  equations,  ,  E  ,  (c)  is  the  self-energy  of  on  electron 

nk 

that  occupies  the  one-electron  Hartrce-Foek  state  r»V  in  an 

(N-l) 

K-electror.  system.  E  (h)  is  the  total  energy  change  in  the 

nk 

remaining  K-l  electron*;  when  the  electron  occupying  the  state  in 

(N-l) 

question  is  removed.  E  (h)  is  the  self-energy  of  a  hole, 

nk 

Fantel'ides  et  al  (  )  have  produced  some  model- 


independent  results  concerning  these  self-energies.  They  find 


that  the  self-energy  of  holes  is  always  positive,  and  therefore 


the  valence  bands  from  the  original  Kartree-Fock  calculation 


always  move  up  on  the  energy  scale  upon  correlation.  These  self¬ 


energies  grow  as  we  approach  the  bottom  of  the  valence  bands,  so 


on  balance  the  valence  bands  arc  narrowed.  The  self-energies  of 


the  lower  electrons  in  the  conduction  band  are  negative,  so  the 


conduction  bands  move  down.  The  models  used  have  only  small 


change  in  the  amount  of  shift  for  the  low-lying  conduction  band, 


and  so  we  use  a  rigid  downward  shift  for  these.  We  end  up 


shifting  the  valence  bands  higher  and  rigidly  dropping  the 


conduction  bands.  This  produces  a  smaller  optical  gap  than  that 


obtained  by  the  Kartrec-Fock  calculations. 


?<> 

We  more  nearly  approach  the  experimental  gap  by  this  narrowing. 

These  self-energies  were  first  calculated  by 

Toyozawa  (  ).  The  theory  of  this  calculation,  the  electronic- 

polaron  method,  was  considerably  further  developed  by  A.  B.  Kunz, 

so  that  it  now  predicts  hole  self-energies  as  well. 

In  the  original  Hartrce-Fock  theory  particles  respond 

only  to  the  average  position  of  the  other  electrons  and  ions. 

This  is  obviously  incorrect,  or  at  least  incomplete?  the 

independent  charge  will  polarize  its  surroundings,  especially  if 

those  surroundings  have  sufficient  time  to  respond.  The  electronic- 

polaron  model  dresses  the  conduction  band  elec trons  and  valence 

band  holes  with  quanta  of  the  polarization  field.  These  quanta  are 

excltons.  In  the  model  we  are  using,  the  excited  states  of  the 

(N) 

crystal  are  simulated  by  a  dispersionlcss  band  of  excltons,  E  (e) 

n  k 

(N-l ) 

and  E  (h)  are  the  interaction  energies  of  a  bare  electron  and 
nk 

hole  with  this  field.  These  energies  are  called  polarization 
energies.  This  model  uses  second-order  perturbation  theory  to 


calculate  the  self-energies. 


V 

These  polarization  energies  are  calculated,  on  the 

"basis  of  a  model  in  which  the  hole  or  electron  is  fixed  in  space. 

This  method,  the  Kott-Littlcton  method,  uses  a  perturbative 
approach  to  find  the  induced  dipole  moments  on  all  ions  of  the 
crystal  due  to  the  localized  charges.  This  calculation  takes  into 
account  the  field  from  the  induced  dipoles  -  it  is  self-consistent. 
Since  it  takes  a  finite  time  for  the  crystal  to  respond  to  such  a 
change,  moving  charges  should  induce  less  polarization  than  the 
model  static  charge.  Therefore,  this  calculation  should  give  an 
upper  bound  for  the  polarization  energies  of  actual,  mobile 
electrons  ana  holes. 

This  electronie-polaron  model  is  basically  a  long- 
range  scheme.  Only  in  the  limit  of  large  distance  can  we  find  these 
changes  by  assuming  that  they  can  be  described  by  dipoles. 

For  short  distances,  the  shape  of  the  wavefunctions  plays  a  major 
role,  and  quadrupole  and  higher  multipole  effects  cannot  be 
ignored.  The  changes  in  the  central  atom  and  its  near  neighbors 


certainly  cannot  be  modeled  by  a  dipole  very  successfully! 


Short-range  correlation  calculates  the  change  of  the 


nearby  orbitals  when  an  electron  Is  added  or  removed.  To  find  these 

corrections,  we  do  simple  atomic  calculations  of  different 

Ionization  states,  supplemented  by  cluster  calc ala tlons  In  the  case 

of  negative  Ions  where  the  added  electron  extends  over  a  significant 

region  of  space.  In  the  limit  of  zero  overlap  between  atoms,  the 

ohort  range  correlation  correction  to  the  hole  energy  is 

ESr(h)  -  E  (4.10) 

scf 

where  <f  is  the  energy  of  the  level  the  electron  was  removed  from 

and  E  is  the  Kartree-Fock  Ionization  energy.  The  ionization 
scf 

energy  is  the  difference  between  the  Kartree-Fock  energy  of  the 
system  with  the  level  occupied  and  the  Kartree-Fock  energy  with 
that  same  level  full.  The  short-range  energies  so  calculated 
are  again  upper  bounds,  since  electrons  and  holes  are  not  as 
localized  as  this  model  potrays  them. 

The  Kartree-Fock  calculation  and  the  following 
correlation  corrections  give  a  set  of  one-electron  energy  bands. 

They  must  now  be  compared  to  experiment  to  test  the  adequacy  of 


our  approximation. 


CADMIUM  SULFID2 


/it  roojn  temperature  and  standard  pressure  CdS 
crystallises  in  the  sine-blende  lattice  with  a  lattice  constant 
of  5*818  angstroms  (10.99  atomic  units)  and  in  the  hexagonal 

ft 

form  with  lattice  constants  of  4,1343  and  6.7234  angstroms  ( 

7,8136  and  12.7054  atonic  units).  The  zinc-blende  form,  the 

subject  of  these  calculations,  is  composed  of  two  interpenetrating 

++  ~  — 

face-centered  cubic  sublattices  occupied  by  Cd  and  S  ions, 

disploced  relative  to  each  other  by  1/4  of  the  diagonal  of  the 

unit  cube.  The  symmetry  properties  of  the  zinc-blencle  lattice 
2 

(  srace  grout*  T  )  have  been  discussed  by  Par. men  ter  (  ) ,  and  in 

d 

the  following  discussions  of  the  band  structure,  and  notation  of 
Kouclreert,  Smoluchowskl  and  Wigner  (  )  will  be  used. 


V«  first  examine  the  energy  states  at  the  Prillouin 

10 

zone  center,  the  point.  Cadmium  ions  contribute  a  filled  (4d) 

shell  to  the  valence  energy  region,  and  sulfur  ions  contribute 

2  6 

a  (3«)  (3p)  configuration.  The  p  and  d  levels  hybridize  to 

some  extent,  since  they  are  separated  by  less  than  nine  evt  still, 

the  valence  bands  are  not  strongly  hybridized.  The  top  valence 

band  is  predominantly  p-llkc  (  ever  75?0  and  the  next  lower  band 

is  predominantly  d-type,  again  over  75#*  We  expect  the  lowest 

conduction  band  to  be  derived  from  the  cadmium  5s  levels,  os  is 

usual  in  compounds  that  have  any  ionic  character.  In  the  zinc- 

blende  structure,  the  crystal  field  splits  the  fivefold  degenerate 

d  levels  into  a  triply  degenerate  state  and  a  doubly  degenerate 

15 

state.  The  threefold  degenerate  p  states  stay  degenerate  under 

12 

this  crystal  field  and  transform  like  .  The  zinc-blende  lattice 

15 

has  no  inversion  symmetry,  so  the  bands  at  the  gamma  point  need 

* 

not  have  a  definite  parity.  The  conduction  band  contains 

15 

both  p  and  d  contributions,  in  fact. 

Several  previous  calculations  have  been  dene.  The  valence 


bands,  especially  the  upper  ones,  ana  the  lower  conduction  bands 


have  been  well-described  by  various  pseudopotentiol  calculations. 


That  is,  they  accord  with  current  experiment.  However,  since 
the  pseudopotential  method  is  basically  a  paramo trlzation  scheme, 
in  which  the  parameters  are  determined  by  experiment,  little 
critical  new  information  can  be  obtained.  If  the  experimental 

evidence  has  been  misunderstood  or  misinterpreted,  the  -  “ 

* 

pseudopotential  will  simply  predict  the  mistakes  or  misunder¬ 
standings  that  it  sprang  fx*om.  A  pseudopotential  calculation 
is  an  aid  .to  undex standing,  but  it  tends  to  not  be  falslflable. 
Such  calculations,  such  as  those  by  Cohen  and  Bergs tresser  (  ) 

also  depend  on  the  validity  of  the  cancellation  theorem,  which 

is  not  exact  and  which  is  much  less  valuable  for  systems  that 

» 

have  localized  states  that  are  closely  comparable  in  energy  with 
the  valence  states  -  here,  the  cadmium  d  states  are  the  ones 
not  easily  described  because  of  their  local  nature. 

A  related  but  more  theoretically  rigorous  technique, 
the  orthogor.alized  plane  wave  method,  has  some  of  the  same  errors. 
In  this  method,  the  valence  and  conduction  states  are  described 


by  plane  waves  orthogonalized  to  the  core  states,  taken  as  constant. 


In  these  previous  CrX  calculations,  as  performed  by  Euwcraa  and 
Stukel  (  ) ,  EuKema  and  Collins  et  al  (  )  and  by  Stukel  et  al  (  ) 

there  have  been  errors  of  up  to  three  ev  in  the  p-like  valence 
and  conduction  bands  in  compounds  made  of  first-row  atoms. 

Also,  the  4d  states  are  misplaced  by  about  a  rydberg  -  probably 
related  to  their  core-like  nature.  The  level  ordering  of  the 
cadmium  4d  and  the  sulfur  s  bands  seems  to  be  reversed  from 
that  measured  by  photoemission  by  about  6  ov.  Altogether  it 
is  easy  to  nee  that  there  is  a  need  for  a  first-principles 
calculation  that  correctly  predicts  the  major  features  of 
the  band  structures  of  these  materials. 

Therefore,  the  first  all-electron,  self-consistent  “ 
nonempirical  band  calculation  using  nonlocal"  exchange  has  been 
performed  on  CdS.  The  emthods  of  calculation  are  as  discussed 
earlier  in  this  work.  It  should  be  mentioned  that  a  selfrconsistent 
band  calculation  has  been  performed  by  Zungcr  and  Freeman  (  ) 

using  a  local  density  approximation  for  the  exchange  potential. 
Their  results,  as  well  as  earlier  ones,  will  be  discussed  in  the 


ensuing  pages. 
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STO  BASIC*  Fl'NCTICKS  FOR  CADSIl’K  A?.*D  LOCAL  OBBIIAL  COEFFICIEKTS 


s  Bails 
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Z 

n 
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41 
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C 

C 

C 

J 

J 

J 
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-.46790 
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O.43795 

0.44531 
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0.33543 

-.00775 

-.01 566 

-.72460 

Baals 

d 

Basis 

Z 

n 

Z 

1J 

2  J 

2  j 

31.29 

3 

20.99 

19.62 

3 

11.48 

10.46 

3 

7 .67 

6.e3 

4 

5.75 

7.00 

4 

2.50 

4.55 

2.90 

2.00 

5 

2.00 
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50 

21 

C 

C 

C 

J 

i 
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0.00337 

-.17167 
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0.10571 
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0.35872 

0.07401 

0.17264 

0.53579 
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-.02815 

1.21004 

-.22539 

0.00925 

-.03467 

0.04950 

-.00394 

-.00954 

0.95921 

32 

42 

C 

C 

J 

J 
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0.30001 

-.23294 
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-.CC499 

O.73322 

O.OC302 
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TALLE  *'/.2 


STO  BASIS  PUHCTIOKS  FOR  SULrlR  AKD  LOCAL  CF.31TAL  COEPFICIERTS 


s  Basis 
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Z  n 

Z 
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OJ  lj 
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3 

1.63  3 

3 

2.32 

1.32 
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-.30626 
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0.02955 
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0.04380  0.00 
-.22859  0.00 
0.54462  0.00 

-.50114  0.00 
0.27367  1.00 


-•46037 

0.34676 


polarizabilities  listed  \  y  Tessron  •'.<*.  ol  (  )  were  then  used  for 

Cd  and  S  along  with  the  onticol  dielectric  constant  of  6.32 

to  calculate  the  polarization  energies  using  the  Kott-Llttleton 

method.  v.fe  obtained  values  of  0.1366  ry  end  0.2534  ry  were 

obtained  os  the  polarization  energies  associated  with  the  S 
++ 

and  Cd  holes  respectively .  I-clng  less  tightly  bound,  the 

sulfur  anions  polarize  to  a  greater  extent  around  a  hole  at  the 
++ 

Cd  cation  site  than  do  c-ndmiur.  ions  around  a  sulfur  hole, 

producing  a  lorgor  polarization  energy  for  the  cation.  We  also 

calculated  short-range  relaxation  cnex'gies  for  the  states  of 

sr  — 

interest  in  tho  valence  region;  these  are  E  (S  )  =  0,07824  ry, 

3P 

sr  —  sr  ++ 

E  (S  )  r=  0. 06302. ry,  and  E  (Cd  )  a  0,116  ry.  These  corrections 
3s  4s 

were  added  to  the  Hamiltonian  matrix  which  was  then  rediagonalized 
to  give  the  correlated  valence  bands.  Conduction  bands  were 

to 

computed  by  shifting  the  Hartree-Fock  conduction  bands  by  -.2584  ry. 


the  polarization  energy  of  a  conduction  band  electron  on  a 
■f+ 

Cd  site. 


Shown  in  figure  4,1  are  the  correlated  energy  bands 


of  CdS.  The  calculated  band  structure  shows  cadmium  sulfide  to  be 


Figure  4,1  Correlated  energy  bands  of  CdS  for  the  normal 


lattice  constant  of  5*^81  angstroms  (  10,994  eu). 


a  direct  hand  rap  semiconductor  with  the  geo  at  the  gamma  point. 

This  conclusion  Is  In  agreement  with  the  previous  band  calculations. 

Three  major  non-overlapping  regions  constitute  the  valence  bands 

in  this  system.  The  S  3s  derived  band  lies  19  eV  below  the 

valence  bond  edge  end  is  about  1,4  eV  wide.  The  next  region 

++ 

is  primarily  derived  from  the  Cd  kd,  and  lies  12  eV  below  the 
top  of  the  valence  bond.  This  band  is  about  one  eV  wide.  The 

t 

«*•  ''t* 

uppermost  valence  band  is  primarily  derived  from  S  3P  levels, 

and  is  about  3-9  cV  wide.  The  lowest  conduction  bond  is  s-llke 

+  + 

and  is  derived  from  the  Cd  5s,  S  3s  and  S  4s  levels. 

The  band  gap  is  found  to  be  direct  and  equal  to  7.1  eV. 

The  optical  value  of  the  gap  ,  E  =  2.55  eV  (  ),  obtained  from 

6 

experiment,  is  in  serious  disagreement,  Thu  correlation  model 
here  used,. the  clectronlc-polaron  model  and  it  limit, 

Mott-Littleton  theory,  is  adapted  to  insulators  and  ignores 
short-range  polarization  effects.  Such  effects  should  be  small 
in  atoms  and  in  systems  where  the  local  orbitals  are  only  slightly 
perturbed  from  the  free  atom  or  ion,  but  in  polar  semiconductors 
such  as  CUS  ,  there  is  no  guarantee  that  the  local  orbitals  are 


Figure  4.2  Density  of’ states  of  CdS 


Figure  ^.3  XFS  spectrum  for  CdS  according  to  Ley  et  al  (  ), 


I  and 
1 


peak.  S 

1 


X  describe  fine  structure  of  the  highest  valence-band 
2 

is  a  shoulder  on  the  high-binding  side  of  this  peak. 


II  represents  the  second  valencc-band  peak. 


that  similar  to  the  free  ions.  In  this  class  of  compounds,  there 


is  substantial  charge  density  in  the  lnternucleor  region. 

It  nos  been  argued  that  such  effects  oct  to  narrow  the  band  gap. 

The  density  of  states  for  CdS  as  calculated  from  the 

original  band  structure  calculation  is  given  in  figure  4.2.  For 

comparison,  figure  4.3  shows  the  corrected  XFS  spectrum  obtained 

* 

by  Ley  et  al  (  ).  The  experimental  evidence  clearly  shows  the 

two  peaks  in  the  density  of  states  of  the  upper  valence  band. 

The  positions  of  these  peaks  are  also  in  reasonable  agreement 
with  experiment.  The  upper  peak  lies  1.6  eV  below  the  top  of 
the  valence  band  by  Ley’s  measurements,  while  ou'  calculation 
gives  a  peak  at  approximately  1.4  e’v.  Experimentally,  the  second 
peak  is  at  4.1  eV,  while  this  calculation  has  a  peak  at  5*^  «V. 

It  also  seems  possible  to  identify  the  shoulder  of  the  upper  peak. 
From  our  band  calculation,  this  shoulder  seems  to  be  at  about 
2.6  eV,  while  Ley"s  measurements  put  it  at  2.1  eV. 

The  position  of  the  d-like  levels  is  also  correctly 
predicted  in  these  calculations.  These  levels,  primarily  formed 
from  cadmium  4d  states,  are  found  to  peak  at9;64  cV  below 
the  top  of  the  valence  band  experimentallly ,  according  to  Ley. 


Figure  4,1/  Joint  density  of  states  of  CdS 


Figure  4.5  Isnaginary  part  of  the  dielectric  constant  of  CdS 

(  this  calculation) 


CdS  IH.  PART  DIELECTRIC  CCRS 


Figure  4.6  Experimentally  derived  imaginary  part  of  dielectric 


constant  in  CdS,  from  Cardona  et  al  (  ). 


Derived  from  ultraviolet  reflection  measurements. 


Cur  calculation  gives  a  value  of  10.°  eV  below  the  tor  of  the 
valence  bend,  in  very  good  agreement  with  experiment.  This  is 

t 

in  sharp  contrast  with  previous  pseudopotcntlal  and  OPW 
calculations,  which  misplaced  this  band  as  much  as  12  cV  below 
its  actual  position.  It  should  be  mentioned  that  Ley's 
photoemission  studies  were  done  on  CdS  in  the  hexagonal  form, 

»  a 

but  since  the  first  two  shells  around  the  central  stom  are 
identical  in  the  hexagonal  and  zinc-blende  structure,  the  atom 
is  exposed  to  a  very  similar  potential,  and  no  substantial  change 
in  the  density  of  states  is  expected. 

In  figure  4.4  the  Joint  density  of  states  as  calculated 
fron  cur  band  structure  is  shown.  In  figure  4.5,  the  Imaginary 
part  of  the  dielectric  function  is  calculated  from  the  some 
theoretical  band  structure.  In  figure  4.6,  Cardona  et  al  extract 
the  imaginary  part  of  the  dielectric  constant  from  ultraviolet 
reflection  measurements  on  zinc-blende  CdS  (  ).  There  is  clear 

qualitative  agreement,  especially  with  the  central  peak  and  the 
two  sub-oeaks  approximately  3  and  5  ®V  above  that  central  peak. 

The  tail  below  the  central  peak  in  energy  is  exaggerated  due  to 
the  overestimate  of  the  band  gap  in  our  calculations. 


To  sum  up:  these  calculations  correctly  predict  the 


cajor  features  of  the  energy  bands  in  CdS,  as  determined  from 
optical  and  other  measurements.  Insofar,  this  first  all-electron, 
self-consistent  non-emplrical  method  has  succeeded.  The  d-like 
levels  arc  correctly  predicted,  os  was  not  the  case  with  earlier 
calculations.  Nothing  in  the  pure  bulk  energy  bands  suggests 
an  explanation  of  the  anomalous  diamagnetism,  ond  it  would  seem 
that  the  effect  is  not  an  intrinsic  one.  It  may  be  hoped  that 
this  better  understanding  of  the  electronic  structure  of  the 
pure  solid  will  aid  in  the  understanding  of  the  anomalous 
diamagnetic  state,  perhaps  by  acting  as  necessary  first  step  in 
the  examination  of  the  properties  of  defects  and  impurities  In 


cadmium  sulfide. 


ZINC  OXIDE 


v* 


At  roo:?i  temperature  and  standard  pressure  ZnO 
« 

crystallizes  in  the  hexagonal  (  wurtzite)  form  with  lattice  constants 

of  3.2^9  angstroms  (6.480  atomic  units)  and  5*19o  angstroms  (9*324 

atomic  units).  The  zinc-blende  form  with  identical  nearest-neighbor 

distance  has  a  lattice  constant  of  4,595  angstroms  (8.634  atomic 

units)  and  is  the  subject  of  these  calculations.  It  is  composed  of 

two  interpenetrating  face-centered  cubic  sublattices  occupied  by 
++■  **•* 

Zn  and  0  ions,  displaced  relative  to  each  other  by  1/4  of  the 
diagonal  of  the  unit  cube. 

We  begin  again  by  examining  the  energy  states  at  the 

10 

Brillouin  zone  center.  Zinc  ions  contribute  a  filled  (3d) 

shell  to  the  valence  energy  region,  and  oxygen  ions  contribute  a 
2  6 

(2s)  (2p)  configuration.  In  this  crystal  structure,  the  threefold 
degenerate  p  states  stay  degenerate  and  transform  like  ' 

15 


The  crysta*  field  splits  the  fivefold  degenerate  dlcvels  into  a 

triply  degenerate  ’  state  and  a  doubly  degenerate  p  state. 

15  12 

The  zinc  3d  levels  and  the  oxygen  2p  levels  lie  quite  close, 
and  this  calculation  shows  shows  all  eight  valence  bands  lying 
within  a  five  eV  region.  The  zinc-blende  lattice  has  no  inversion 
symmetry,  so  the  bands  at  the  gamma  point  need  not  have  any 
definite  parity.  This  calculation  shows  that  this  upper  valence 
region  exhibits  strong  p-d  hybridization.  The  band  lying 
beneath  this  upper  valence  region  lies  nearly  24  eV*  below  the 
top  of  the  valence  band,  and  is  primarily  derived  from  oxygen 
2s  states.. 

In  general,  in  ionic  compounds  such  as  ZnO  we  would  expect 
the  lowest  conduction  band  to  be  primarily  formed  from  the 
zinc  4s  level.  The  actual  calculation  shows  that  although  the 
zinc  4s  states  play  a  major  role,  oxygen  2s  and  3s  states  actually 
play  a  larger  role  in  this  first  conduction  band.  The  next, 
conduc t ion.  band, ■■  tx*i ply  degenerate  at  the  ".one  center,  are  formed 
almost  entirely  from  zinc  4p  states  .  ‘This  is  suggested  merely 
by  comparing  zinc  with  its  successor  in  the  periodic  table,  gallium. 


Several  previous  calculations  liave  been  performed  on 


zinc  oxide.  Attempts  using  local  pseudopotential  theory,  such  os 
those  by  Rossler  (  )  and  Bloom  and  Ortenburger  (  )  have  not 

been  entirely  successful.  Since  the  pseudoootential  theory  depends 
upon  a  cancellation  cf  the  strong  core  part  of  the  potential  by 
the  usual  requirement  that  the  valence  electrons  be  orthogonal 
to  the  core  electrons,  first-row  elements  would  seem  to  be 
unsuitable  for  this  approach.  The  point  is  that  valence  p-states 
for  first  row  elements  are  not  required  to  be  orthogonal  to  any 
p-core  states.  Investigators  have  attempted  to  alleviate  this 
problem  by  empirical  adjustments  to  the  pseudopotentlols,  but 
it  hasn’t  worked  well.  When  nonlocal  pseudopot-nntials  were  used, 
as  in  the  calculations  of  Phelikowsky  (  ) ,  a  better  agreement 

with  experiment . is  achieved. ‘Even  so,  much  of  the  value  of  the 
pseudopotential  approach  is  lost.  There  is  not  the  ssme  confidence 
that  the  pseudopotential  will  retain  its  preditetive  powers  in 
different  compounds,  because  the  physical  rationale  is  weakened. 
The  modified  and  nonlocal  calculations  still  suffer  from  all  the 
original  weakness  of  pseudopotentials:  reliance  on  experiment. 


STO  FAS  IS  FUNCTION’S  FOR  ZXKC 


AND  LOCAL  ORBITAL  COEFFICIENTS 


s  'Basis  p  Eaais  d  Basis 
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TABLE  CP 

STC  BASIS  FUNCTIONS  FOP.  OXYGEN  AND  LOCAL  ORBITAL  COEFFICIENTS 


a  Basis  p  Basis 
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z 
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-.17832 
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1.45654  -.68094 
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0.03913 
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0.45277  -.06106 
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0.14249 

0.71 512 

-.00728  0.01303 
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-.00830 

0.49634 

-.43478 

2 

0.04136 

-.32320 

0.22209 

Therefore,  we  perform  an  all-electron,  ftnlf-conssistent 
noneraplricnl  band  calculation  on  ZnO.  The  methods  of  cumulation 
ere  the  same  as'  used  previously  on  CdS. 

The  calculation  was  besun  by  colculeting  local  orbitals 
for  zinc  oxide  with  the  normal  lattice  constant.  The  zinc,  and 
oxygen  basis  sets  of  Ragus  et  al  (  )  were  used  and  were  modified 

for  this  calculation.  Two  additional  diffuse  STOs  were  added  to  the 
s  basis  on  zinc.  The  results  of  the  local  orbitals  calculation 
are  listed  in  tables  5*1  and  5.2  for  zinc  and  oxygen  ions, 
respectively, 

Hartree-Fock  bands  were  calculated  for  20  k  points 

in  the  irreducible  wedge  of  the  Brillouin  zone.  The  usual- 

exaggerated  band  gap  appears?  we  proceed  to  apply  the  correlation 

methods  previously  discussed. 

Polarizabilities  listed  by  Tessman  et  al  (  )  were*  then 

+*f 

used  for  Zn  end  0  along  with  the  optical  dielectric  constant 

of  4.036  to  calculate  the  polarization  energies  using  the  Kott- 

Littleton  method.  We  obtained  values  of  0.2172  ry  and  0.2364  ry 

++ 

as  the  polarization  energies  associated  with  the  0  and  Zn 


holes  respectively.  We  also  calculated  short-range  relaxation 


Figure  5.1  Correlated  energy  bands  of  Zno  for  the  normal 


lattice  constant  of  4,595  angstroms  (  8.634  au) 


energies  for  the  state  of  Interest  in  the  valence  region i  this 
sr  ++ 

As  E  (Zn  }  =  0.31^2  ry.  These  corrections  were  added  to  the 

3d 

Hamilton  Aon  matrix  which  was  then  redlagonallzed  to  give  v.*e 

correlated  valence  bands.  Conduction  bands  were  computed  by  shifting 

the  Kartree'Fock  conduction  bonds  by  -,2864  ry,  the  polarization 

++ 

energy  of  a  conduction  bond  electron  on  a  Zn  site. 

Shown  in  figure  5,1  are  the  correlated  energy  bonds 

of  ZnO.  The  calculated  bond  structure  shows  zinc  oxide  to  be 

a  direct  bond  gap  semiconductor  with  the  gap  at  the  gamma  point. 

This  conclusion  is  in  agreement  with  previous  band  calculations. 

Two  mojor  nor.-overlapplng  regions  constitute  the  valence  bands  in 

this  system.  The  highest  region  is  about  4  cv  in  width,  and  is 

divided  into  two  subregion^' with  almost  no* overlap.  The  higher 

end  wider  of  these  subregions  is  ,  approximately  3  eV  in  width, 

—  ++ 

is  composed  of  0  2p  and  Zn  3d  levels,  while  the  lower  subband 

++ 

is  almost  dispersionless  and  is  of  almost  pure  Zn  3d  character. 

t 

The  lowest  valence  band  is  about  2.3  eV  in  width  and  is  composed 
of  0  2p  and  Zn  3d  levels. 

The  band  gap  is  found  to  be  direct  and  equal  to  10.86  eV. 
Since  the  gap  is  experimentally  found  to  te  3.3  eV,  it  must  be 


Figure  5*2  Density  of  states  of  ZnO 
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Flsure  5.3  Corrected  XFS  spectrum  for  2n0  according  to  Ley  et  el(  ). 
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that  the  conduction  band  structure  here  predicted  5s  seriously 
in  error.  Previously  mentioned  potential  problems  in  our  correlation 
model  nay  explain  this  discrepancy,  but  it  seems  most  liV.ely  that 
there  is  some  flaw  in  the  estimate  of  the  polarization  energies'. 

The  density  of  states  for  Zn 0  as  calculated  from  our 
band  structure  is  given  in  figure  5,2.  For  comparison,  figure  5.3 
shows  the  corrected  experimental  density  of  states,  derived 
from  x-ray  photoemission  experiments  conducted  by  Ley  et  al  (  ), 

The  d-like  levels  have  been  subtracted  out  from  the  experimental 
data,  since  they  dominate  the  spectra.  The  double  peak  in  the 
density  of  states  of  the  uppermost  band  is  clearly  shown  ir. 
experiment  and  in  our  calculation.  The  d-like  levels  in  our 
i  calculation  have  some  structure,  are  not  smoothed  into  a  single 
.peak  as  in  the  experiment.  The  experimental  peak  is  centered 

C  . 

‘.8,81  eV  below  the  top  of  the  valence  band,  according  to  Ley. 

Our  d-complex  is  situated  approximately  5  eV  below  the  top, 
but  the  qualitative  picture,  that  of  a  double-peaked  valence 
band  with  a  very  sharp  d-band  .about  3  eV  below,  corresponds 


closely  with  experiment. 


Our  calculations  also  (jive  «  uand.  derived  from  0  2a 


states  (not  shown  on  figure  5.2)  centered  around  eV  below 
the  top  of  the  valence  band  -  this  corresponds  to  the 
peak  in  figure  5,3 

In  addition,  we  have  calculated  the  joint  density  of 
states  (  figure  5«M  ««d  the  imaginary  part  of  the  dielectric 
constant  (  figure  5*5)  for  zinc  oxide. 


Figure  5*//  Joint  Density  of  States  of  ZnO 
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Figure  5. T Imaginary  part  of  the  dielectric  constant  of  ZnO 


CONCLUSIONS 


We’ve  used  ab  Initio  Hartree-Fock  theory,  as  well  as 

relaxation  and  polarization  correlation  corrections  to  that  theory, 

to  calculate  the  electronic  band  structures  of  cadmium  sulfide  and 

zinc  oxide.  In  the  calculation  of  cadmium  sulfide,  hybridization 
++  ~  —  —  - 

of  the  Cd  5s,  S  3s  and  S  bs  is  a  major  factor  in  the  lowest 
conduction  band.  The  calculation  shows  that  CdS  is  a  direct  gap 
semiconductor  with  a  gap  of  7*1  eV  at  the  gamma  point.  This  tends 
to  show  that  the  mechanism  proposed  by  Abrikosov  to  explain  the 
anomalous  diamagnetism  originally  seen  in  CuCl  (  and  discredited 
there  by  Weldman  (  )  )  is  not  applicable  to  CdS,  since  it  requires 

a  small  indirect  gap. 

The  calculated  one-electron  energy  bands  are  compared 
with  the  published  optical  data.  The  valence  bands  are  in  excellent 


agreement  with  photoenission  and  reflection  data,  and,  in  particular 


the  position  of  the  d-llke  band  is  correctly  predicted,  in  contrast 
with  previous  calculations. 

The  calculation  for  zinc  oxide  shows  that  there  is 
very  significant  nixing  in  the  uoper  valence  bands  between  the 
0  2p  and  Zn  3d  levels.  It  seems  clear  that  previous  pseudo potential 
calculations  that  could  not  correctly  take  into  account  the  zinc 
d  levels  were  incanable  of  explaing  the  valence  bands.  The  calculation 
gives  a  band  gap  of  10.86  ev,  and  predicts  ZnO  to  be  a  direct  gap 
semiconductor  with  the  gap  at  the  gamma  point.  In  the  case  of  both 
CdS  and  ZnO,  it  should  be  recalled  that  we  are  essentially  Just  solving 
Dyson's  equation,  and  that  higher  order  diagrams  are  quite  capable 
of  accounting  for  the  difference  between  the  (exaggerated)  gaps 
from  our  calculations  and  the  measured  experimental  gaps.  The  fact 
remains  that  these  calculations  predict  the  valence  structure? 
the  one  electron  bands  for  zinc  oxide  also  accord  well  with  published 
optical  experiments.  Again,  the  position  of  the  d-like  bands  is 
correctly  predicted. 

The  anomalous  diamagnetism  observed  in  CdS  remains 
unexplained.  It  seems  to  be  an  extrinsic  phenomenon  . 


The  observation  of  ferromagnetism  upon  an  applied 


field,  as  well  as  the  anomalous  diamagnetism,  is  intriguing. 

It  suggests  that  the  pairing  mechnaisn  may  favor  a  state  with  a 
spin  of  one,  as  opposed  to  the  Cooper  pair,  with  a  net  spin  of  zero. 
Each  pair  would  have  an  intrinsic  magnetic  moment.  This  possibility 
has  been  discussed  for  classical  phonon-mediated  superconductivity 
by  ?.  W.  Anderson  and  P.  Korel(  )  but  seems  not  to  have  been 
observed,  except  possibly  at  extremely  low  temperatures  (  ). 

An  exciton-medlated  pairing  might  well  favor  the  1*1 

or  higher  state.  The  superfluid  state  in  He  is  suggestive. 

3 

As  yet  the  details  of  the  interactions  are  not  know, 
although  it  does  seen  that  chemical  impurities  play  a  role.  It  may 
still  be  possible  to  test  some  of  these  ideas  by  a  phenomenological 
theory  like  the  Glnzburg-Landau  theory.  One  possibility  seems 
interesting*  in  the  Ginzburg  -Landau  theory,  two  characteristic 
lengths  appear,  the  coherence  length  and  the  penetration  length. 

In  such  a  model  revised  for  a  p-wave  system,  a  third  length  arises, 
the  characteristic  length  for  the  change  in  spin  direction. 

The  two  lengths  of  Ginzburg-Landau  theory  allow  the 


existence  of  a  distinctive  surface,  which  may  have  a  lower  energy 


than  the  bulk  state,  which  leads  to  the  interesting  and 


technologically  useful  type  IX  superconductors.  This  still  exists 
with  three  characteristic  lengths,  but  another  effect,  becomes  possible* 

A  second  surface  layer  is  introduced;  we  may  compare  this  to  the 
earth,  with  crust,  mantle  and  core.  If  the  mantle  is  energetically 
favored,  bubbles  of  a  certain  sire  would  be  energetically  favored. 

These  superconducting  domains  might  explain  the  very  high  but  finite 
conductivity  seen  in  the  anomalous  states  of  CuCl  and  CdS,  and 
alignment  of  these  domains  might  explain  the  fexromagnetlsm  seen 
in  CdS  under  high  applied  field.  At  the  moment,  this  is  all  speculation. 
We  await  further  experiments. 
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ABSTRACT 

Energy  bands  and  ionization  energies  for  electrons  in  A1P,  GaP,  and 
CaAs  are  obtained  using  a  new  method  proposed  by  Kunz,  et  al. ^  The  valence 
bands  and  ionization  energies  obtained  are  found  to  be  in  good  agreement 
with  experiment.  The  conduction  bands  vary  substantially  from  experiment. 
Finally,  a  method  for  obtaining  better  conduction  bands  is  proposed. 


1.  INTRODUCTION 


A  new  ab  initio  method  for  obtaining  correlated  energy  bands  was 
recently  proposed  by  Kunz,  et  al.^  The  need  for  a  new  method  becomes 
apparent  when  the  inherent  problems  associated  with  current  ab  initio 
methods  are  considered.  Existing  methods  can  be  broken  into  two  classes, 
those  which  use  the  nonlocal  Hartrec-Fock  (H-F)  exchange  potential  and 
those  which  use  some  local  density  approximation  to  the  exchange  or 
exchange-correlation  potential.  The  former  class  of  methods  requires 
calculations  which  arc  both  lengthy  and  complicated.  Once  such  a  H-F 
calculation  has  been  performed,  there  is  still  a  need  to  obtain  correlation 
corrections.  These  corrections  are  poorly  understood  for  covalent 
semiconductors  like  the  III-V  compounds.  On  the  other  hand,  local  density 
approximation  calculations  give  results  which  are  highly  dependent  cn  the 
choice  of  the  exchange-correlation  potential.  No  single  potential  has  been 
found  to  give  good  results  for  all  classes  of  compounds.  For  example, 
Herman,  et  al.2,  found  that  in  using  both  the  Slater^  and  the  Kohn-Shan*1 
exchange  potentials  for  several  semiconductors,  the  Slater  potential  agreed 
more  closely  with  experiment  for  II-VX  semiconductors,  while  Kohn-Sham 
exchange  gave  better  results  for  the  III-V  compounds.  Even  in  cases  where 
a  given  potential  yields  the  correct  valence  and  conduction  band 
structures,  such  calculations  have  not  always  been  able  to  place  the  core 
levels  in  their  correct  locations  relative  to  the  valence  bands.  Finally, 
local  density  calculations  have  not  been  successful  ir.  obtaining  the 
ionization  energies  for  electrons  ir.  most  compounds.  We  believe  that  many 
of  these  difficulties  are  a  result  of  using  a  Hartree  potential  which 
includes  a  self-repulsion  tern,  on  the  assumption  that  the  potential  from 
one  electron  in  an  infinite  crystal  will  be  negligible.  As  has  been  shown 
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by  Kunz,  ct  ol.  this  assumption  is  not  correct  in  many  cases  and  is  a 
source  of  significant  error  for  insulators. 

Tnc  new  ncthod,  which  will  be  referred  to  as  the  Hartrce-plus  method, 
uses  the  correct  Hartroe  potential,  _i.e.  with  the  self-repulsion  removed, 
along  with  a  local  exchange-correlation  potential.  The  Hartrce-plus  method 

ft 

was  applied  by  Kunz,  et  al.^t  to  the  solid  rare  gases  and  N’aCl  yielding 
bands  which  are  in  good  agreement  with  both  experiment  and  previous  H-F 
plus  correlation  calculations.  Tne  III-V  semiconductors  provide  a  good 
test  of  the  range  of  applicability  for  the  Hartrcc-plus  method  since  they 
arc  wide  band  semiconductors,  r.s  opposed  to  the  rare  gases  ar.d  HaCl  which 
are  narrow  band  insulators. 

In  Sec. 2,  H-F  theory  is  developed  and  correlation  i3  discussed.  The 
Hartrce-plu3  theory  and  details  of  the  calculations  are  presented  in  Sec. 3 
and  See Ji  respectively.  In  Sec. 5,  the  results  of  Hartree-plus  bar.c 
calculations  for  A1P,  GaP,  ar.d  GaAs  arc  presented  and  compared  to 
experiment  and  previous  theoretical  calculations. 

2.  BASIC  THEORY 


Tne  system  of  interest  contains  n  electrons  and  N  nuclei.  Using  the 
Born-Oppenhcimer  approximation^,  ar.d  neglecting  relativistic  effects,  the 
Hamiltonian  is: 


H-  -  £  jw  i!  w 


1 1 1 


wnere: 


// 


•  V'  in-Ril 

_  -  3. 

?y "  r*-$i 

The  energy  is  in  Rydbergs,  e2=2,  uppercase  letters  refer  to  nuclei, 
lowercase  letters  denote  electrons,  and  the  prime  on  the  sum  indicates  that 
the  self-repulsion  term,  i=j,  is  excluded.  The  problen  now  is  to  solve  the 
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Schrodinger  equation: 


(X,  •••**) 


where  is  the  exact  n-electron  wavefunction  for  the  ith  excited  state, 
and  the  are  taken  to  include  both  space  and  spin  coordinates. 

To  simplify  this  problem,  we  approximate  the  exact  wavcfur.ction 
with  a  single  Slater  determinant  formed  from  a  set  of  onc-electron 


orbitals,  (^(xj)),  i=1,***,n  and  j=1,**‘,n: 

-  fyj)  &(*,)]} 

The  expectation  value  of  H'"',  for  the  trial  function  will  be  a 

rigorous  upper  bound  to  the  exact  ground  state  energy  if  the  resulting 
energy  El*'t>  is  stable  against  variation  in  the  j^.  Performing  such  a 
variation,  subject  to  the  constraint  that  the  r>x  be  orthonornal,  yields  an 
equation  for  the  ^ : 

F#  (x,)  “  f  A,y  <P,CK,) 

where  F  is  the  Fock  operator: 

F=  f,  -pa,  x‘y>p(2/  ol  ( 3 ) 

pcx.xi)-  zcc  fya.wfaj) 

In  eq.(3),  P(2,1)  is  an  exchange  operator,  the  nay  be  used  as 

operators,  znd  is  the  Fock-Dirac  density  matrix.  The  second  term 

is  known  as  the  direct  or  Hartree  terra,  and  the  last  term  is  the  exchange 
term.  Note  that  the  prime  is  no  longer  necessary  since  the  Hartree 
self-energy  is  .exactly  cancelled  by  the  self-exchange.  Diagonalizing  the 
.  matrix Aij  gives  the  standard  Hartree-Fock  equation: 

Fti-a,)  =  €?  fra,)  w) 


The  eigenvalues  of  eq.(4)  are  given  meaning  by  Koopmans'  theorem^. 
Labelling  occupied  states  with  1=1,  •••,n  and  virtual  states  with  a=n+1,‘“, 
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the  eigenvalues  arc: 

«=;= 

<a  =  er0”-  er° 


where  E^* is  the  energy  of  the  single  Slater  determinant  approximation  to 
the  total  m-electron  wavefur.ction  obtained  by  adding  or  removing  the  «th 
orbital  to  the  n  orbitals  used  to  form  the  H-F  ground  state.  Thus,  and 
£,  are  the  energies  that  an  electron  would  have  if  its  instantaneous 
location  were  independent  of  the  instantaneous  locations  of  the  other 
electrons,  and  if  no  relaxation  were  allowed.  To  improve  on  the  H-F 
eigenvalues  while  retaining  the  independent  orbitals,  one  defines 
correlated  energies  by  replacing  the  approximate  energies  in  eq.(5)  with 
the  exact  energies: 


W 


Cm  "  Ctt  Cm 


The  energies  G;  and  are  oua£iparticle  excitation  energies.  The 
quasiparticles  associated  with  these  energies  are  referred  to  as  holes  and 
conduction  electrons  respectively. 

There  are  two  basic  approaches  to  finding  the  quasiparticle  energies. 
The  better  understood  of  the  two  is  the  H-F  plus  correlation  method,  which 
is  essentially  a  perturbative  approach.  In  this  method,  eq.(*0  is  solved 
first,  and  then  the  correlation  effects  arc  added  in  as  a  small  number  of 
corrections.  This  approach  to  obtaining  correlated  bands  was  reviewed  and 
developed  in  a  formal  way  by  Pantelides,  et  al.7  The  correlation 
.  corrections  are  fairly  well  understood  for  both  metals  and  narrow-band 
wide-gap  materials.  Unfortunately  the  III-V  compounds  do  not  fall  into 
either  of  these  classes.  The  other  approach  is  to  use  an  equation, 
analogous  to  the  H-F  equation,  which  somehow  incorporates  correlation  into 
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the  potential.  This  is  the  approach  taken  by  the  Har tree-plus  method. 

3.  THE  HAM'HEE-PLUS  METHOD 

In  the  past,  band  calculations  using  local  exchange-correlation 
potentials  have  solved  an  equation  of  the  form: 

Os.r.  ba,)  -  e{  4>}  (X,)  (?) 

where: 

0$f,  c/7>  p  (•***>,)  +  Vxc  Iptr)  J 

is  a  single  particle  operator  which  replaces  the  nonlocal  exchange  in  the 
Fock  operator  with  a  function,  Vxc  ,  of  the  local  density,^?).  Unlike 
the  Fock  operator,  cq.(7)  no  longer  has  an  exact  cancellation  of  the 
self-repulsion  by  a  self-exchange.  It  has  generally  been  assumed  that  the 
self-repulsion  is  negligible  since  it  come3  from  one  electron  in  an 
infinite  crystal.  This  is  in  fact  true  for  an  electron  in  a  totally 
delocalized  Bloch  state.  However,  what  i3  actually  desired  is  the  energy 
required  to  add  one  quasiparticle  to  the  ground  state  system.  For  all 
insulators,  and  many  semiconductors,  these  quasiparticle  orbitals  are 
local.  In  this  case  it  is  reasonable  to  work  in  the  local  Heitlcr-London 
representation.  For  such  local  orbitals,  the  self-energy  is  not  in  fact 
small  and  should  be  removed  explicitly. 

In  the  local  representation,  let  #*(r)  be  the  i^*  local  orbital  at 
site  M.  The  corrected  potential  then  is: 

-f,  +■  +  VK  V) 

where: 

VJ  C?,)  -  y,*  ptXzXl  )  —  fsfe 

is  the  nonlocal  Hartree  potential.  Although  the  Hartree  potential  is 
nonlocal,  the  nonlocality  is  easy  to  deal  with  by  using  a  single  Hartree 
potential  for  all  orbitals  of  the  form: 

<4,,,!  (fj 


where: 


Zf  ( S.S. )  S  2  /  W,ulZ)l '  l$«(ii)lx  ln.-%l  -'Jh  <** 

Here  the  bra-ket  notation  is  used  to  emphasize  the  operator  nature  of  the 
self-repulsion  term. 

For  Vxc  we  have  used  the  local  exchange  potential  of  Kohn-Shan1*  and 
the  dielectric  screening  function  of  Robinson,  Bassani,  Knox,  and 
Schrieffcr(RBKS)8.  Thus: 

Y/c  LP<»1  =  Vfir)  Fr'“<«)  (J  o) 

where: 

v<um  =  -‘t[(^?)/cc^T'3 

FM = t~  "5  *  V«  0^ 

and: 

Although  it  might  seem  inconsistent  not  to  remove  the  self-exchange  in 
eq.(ll),  one  can  argue  that  the  self-exchange  will  only  be  a  significant 
fraction  of  the  total  exchange  potential  in  regions  of  low  density  where 
F(s<.)  approaches  zero, 
h.  THE  CALCULATIONS 

To  calculate  the  Hartrce-plus  bands  shown  in  Figs.  1-3,  we  first 
performed  an  LCAO  calculation  using  H-F  atomic  wavefunctions  obtained  by 
Bagus,  Gilbert,  and  Roothaan^.  in  this  calculation  Nesbet's  symmetry  and 
equivalence  restrictions1®  were  applied.  Also,  the  crystal  potential  was 
approximated  by  the  first  term  of  its  expansion  in  Yjn,  i.e.  1=0.  The 
next  tern  in  the  potential  would  be  1=3  and  should  not  have  a  large  effect 
on  the  s,  p,  and  d  orbitals,  which  are  all  that  are  of  interest  here.  The 
potential  and  matrix  elements  are  evaluated  using  10  shells  about  the 
central  site.  L'e  represent  the-  remainder  of  the  crystal  by  a  residual 


Kadelung  potential.  To  obtain  this  potential,  we  used  effective  charges 
given  by  Phillips1 1,  This  LCAO  calculation  can  not  be  expected  to  give 
good  results  for  the  conduction  levels  due  to  the  use  of  Hesbet's 
restrictions10.  The  self-energies  of  the  atonic  orbitals  have  been  removed 
for  both  the  occupied  and  the  unoccupied  parts  of  the  orbitals.  This  means 
that  the  conduction  energies  have  been  obtained  as  though  the  conduction 
electrons  felt  the  potential  from  n-1  electrons,  instead  of  all  n  valence 
electrons.  To  remedy  this  the  self-energies  are  explicitly  reinserted  for 
the  virtual  states.  At  this  point  the  Hartrce-plus  matrix  is  no  longer 
diagonal.  However  the  coupling  between  the  occupied  and  virtual  states  is 
very  weak  and  a  second  diagonalization  produces  new  valence  bands  which, 
for  GaP,  differ  in  energy  from  the  previous  bands  by  at  most  HJ.  Taking 
all  this  into  consideration,  we  believe  that  the  final  Hartree-plus 
operator  so  obtained  is  consistent  with  the  configuration  used.  Finally  we 
form  a  Hartree-plus  matrix  using  atonic  orbitals  for  the  core  states  and  89 
planewaves  for  the  conduction  and  valence  states,  and  diagonalize  it  to 
obtain  the  final  band  structure.  For  A1P,  we  repeated  this  calculation 
using  27,  51,  and  65  planewaves. 

5.  RESULTS 

In  all  the  band  calculations  performed  here,  energy  levels  were 
obtained  at  21  nonequivalent  points  in  the  Brillouin  zone  for  an  f.c.c. 
lattice.  The  bands  are  drawn  along  axes  of  high  symmetry  which  connect  the 
points  of  high  symmetry  in  the  Brillouin  ;:one.  The  symmetry  labels  used 
are  consistent  with  Parnenter1^.  In  all  cases  the  origin  has  been  chosen 
to  be  the  cation  site.  In  addition,  the  orbitals  and  their  energies  are 
superscripted  to  indicate  whether  they  are  valence  (v)  or  conduction  (c) 


states. 
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The  valence  band  structures  for  A1P,  obtained  using;  27.  51,  65,  and  89 
planewaves^  are  summarized  in  Table  1.  Mo  experimental  or  theoretical 
values  were  available  for  comparison.  The  bands  obtained  with  09 
planewaves  are  shown  in  Fig.l.  By  comparing  the  results  for  several 
numbers  of  planewaves,  we  can  estimate  the  size  of  the  errors  due  to 
incomplete  convergence.  For  the  top  two  p-derived  bands,  our  results 
appear  to  be  converged  to  within  0.1  cV.  The  third  p-like  band  is 
converged  to  within  about  0.3  eV,  and  the  separation  of  the  s~band  and  the 
p-bands  is  converged  to  within  roughly  O.H  cV.  On  the  otherhand,  the  basic 
structure  is  stable,  as  is  the  width  of  the  s-band.  We  assume  that  this 
convergence  information  is  also  valid  for  our  other  calculations. 

In  Tables  2a  and  2b,  the  valence  band  structures  for  GaP  and  GaAs 
shown  in  Figs. 2  and  3  are  compared  to  experiments  and  previous  theoretical 
calculations.  It  is  particularly  interesting  to  compare  the  Hartree-plus 
results  with  the  two  first  principles  OPW  calculations.  For  both  GaP  and 
GaAs,  the  Hartree-plus  method  gives  valence  band  widths  which  arc  in  better 
agreement  with  experiment  than  those  obtained  in  the  first  principles  OPW 
calculations,  and  are  comparable  ir.  quality  to  the  empirically  adjusted 
calculations,  especially  when  the  lack  of  convergence  is  taken  into 
consideration.  All  of  the  other  calculations  give  better  results  than  this 
work  for  the  separation  of  the  s-band  from  the  p-bands.  This  error  ir.  the 
Hartree-plus  results,  1.0  to  2.0  eV,  is  probably  due  to  incomplete 
convergence  and  relaxation  effects  that  have  not  been  included.  None  of 
the  other  calculations  listed  give  the  position  of  the  d-bands  relative  to 
the  top  of  the  valence  bands.  For  the  Hartree-plus  results,  we  believe 
that  the  errors  in  the  d-bands,  which  are  less  than  15/5,  are  a  result  of 
relaxation  effects,  which  would  be  expected  to  shift  the  d-bands  up 
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relative  to  the  top  of  the  valence  bands. 

In  Table  2c,  the  ionization  energies  for  the  cation  d-levcls  of  GaP 
and  GaAs  are  compared  to  experimental  values.  This  information  is  not 
given  by  any  of  the  previous  band  calculations  mentioned.  The  differences 
between  theory  and  experiment  are  easily  accounted  for  by  relaxation  and 
surface  polarization  effects. 

In  Table  3,  the  lowest  transition  energies  are  compared  to  experiment. 
Clearly  the  Hartree-plus  conduction  bands  are  in  very  poor  agreement  with 
experiment.  To  test  for  errors  in  the  treatment  of  the  self-energies,  we 
performed  a  calculation  for  A1P  with  the  self-energies  set  equal  to  zero. 
Tnis  calculation  yielded  conduction  bands  which  differ  from  those  with 

nonzero  self-energy  by  no  more  than  0.2  eV.  We  therefore  conclude  that  the 

errors  in  the  conduction  bands  are  due  to  the  choice  of  the 

exchange-correlation  potential ,  VXc-  It  is  not  really  surprising  that  our 
choice  of  exchange-correlation  potential  gives  good  results  for  the  valence 
bands  but  not  for  the  conduction  bands.  As  was  shown  by  Pantelides,  et 
al.?,  for  both  insulators  and  wide  gap  semiconductors,  correlation  effects 
in  valence  bands  are  a  result  of  virtual  scattering  of  holes,  while 
correlation  of  the  conduction  bands  cones  primarily  from  virtual  scattering 
of  electrons.  Since  the  mechanisms  for  correlation  in  valence  and 
conduction  bands  are  different,  it  is  unreasonable  to  expect  one 

exchange-correlation  potential  to  be  valid  for  both  cases.  We  suggest  that 
in  future  Hartree-plus  calculations,  it  should  be  possible  to  obtain  better 
results  by  using  two  separate  potentials  for  the  valence  and  conduction 
states. 

6.  CONCLUSION 

Based  upon  the  results  given  by  Kunz,  et  al. ^ ,  for  the  rare  gases  and 
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the  results  for  III-V  semiconductors  reported  here,  we  conclude  that  the 
Hartree-plus  method  should  give  good  results  for  the  valence  bands  of  both 
insulators  and  semiconductors  when  used  with  the  R3KS  screened  Kohn-Shaa 
exchange.  This  method  not  only  gives  good  results  for  the  valence  band 
widths,  it  also  gives  reasonable  values  for  the  core  d  levels  and  the 
ionization  potentials.  Based  upon  the  calculations  given,  we  suggest  that 
it  might  be  possible  to  find  a  single  exchange-correlation  potential  that 
could  be  used  to  obtain  good  conduction  bands  for  insulators  and 
semiconductors  in  a  two-potential  band  calculation.  Such  a  calculation 
would  use  one  potential  for  valence  states  and  a  second  potential  for 
conduction  states.  This  type  of  calculation  not  only  would  be  expected  to 
give  better  Hartree-plus  band  structures,  but  would  also  be  more  consistent 
with  what  we  know  about  correlation  in  semiconductors  and  insulators. 
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CHAPTER  1 


INTRODUCTION 


The  electronic  structure  of  surfaces  is  important  in  many 
scientific  and  technological  fields,  such  as  corrosion,  semiconductor 
devices,  and  catalysis.  However,  the  understanding  of  fundamental  solid 
surface  properties  has  lagged  behind  the  progress  in  understanding  of 
bulk  properties.  This  situation  has  been  a  result  of  the  lack  of 
experimental  methods  for  the  quantitative  characterization  of  the 
surface  structure,  and  the  theoretical  complexities  resulting  from  both 
the  loss  of  symmetry  at  the  surface  and  any  possible  reconstruction 
occuring  at  the  surface.  Recently  developed  experimental  techniques 
such  as  electron  energy  loss  spectroscopy  (ELS)  (  ),  two  photon 
picosecond  spectroscopy  (  ),  surface  extended  x-ray  absorption  fine 
structure  (EXAFS)  (  ),  and  extended  appearance  potential  fine 
structure  (EAPFS)  (  )  have  made  accessible  the  experimental  study  of 
the  structural  and  electronic  properties  of  real  surfaces.  The 
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development  of  powerful  computational  methods  has  made  possible  the 
theoretical  study  of  these  systems. 


1.1  THE  INITIAL  ADSORPTION  OF  OXYGEN  ONTO  THE  ALUMINUM  (100)  SURFACE 


Aluminum  surfaces  oxidize  readily  upon  exposure  to  molecular 
oxygen.  Several  experimental  studies  have  been  performed  on  the  initial 
oxidation  of  aluminum  surfaces  (  ).  Ultraviolet  photoemission 
spectroscopy  (UPS)  (  )  and  Auger  spectroscopy  (  )  studies  have  been 
reported  for  many  crystal  orientations  of  aluminum,  including  the  (100) 
surface.  LEND  (  ),  EXAFS  (  ),  and  EAPFS  (  )  studies  have  been  reported. 

Ultraviolet  photoemission  experiments  (  )  have  reported  a  valence 
band  resonance  at  1.5  eV  below  the  Fermi  level  for  low  oxygen  exposure. 
This  resonance  is  attributed  to  the  oxygen  2p  band.  The  shift  of  the 
aluminum  2p  core  level  has  also  been  measured  for  both  low  and  high 
oxygen  exposures.  This  shift  toward  lower  energy  is  due  to  the  dipole 
moment  arising  from  the  charge  transfer  from  the  aluminum  substrate  to 
the  adsorbed  oxygen.  Work  function  measurements  have  also  been  reported 
(  ),  with  the  work  function  found  to  decrease  with  increasing  oxygen 
exposure.  This  decrease  has  been  interpreted  as  being  a  result  of 
incorporation  of  the  electronegative  oxygen  adatoms  beneath  the  aluminum 


surface. 
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LEED,  EXAFS,  and  EAPFS  techniques  are  able  to  probe  the  surface 
geometry  of  the  aluminum  plus  oxygen  adsorbate  system,  and  thus  these 
studies  are  of  special  interest  as  they  will  help  to  differentiate 
between  various  theoretical  models.  For  low  oxygen  coverage,  the 
reported  Al-0  bond  length  is  1.98  8  (  ),  in  agreement  with  previous 
theoretical  calculations  (  ).  For  higher  exposures,  a  value  of  1.88  A 
is  reported  (  ),  which  correspondsto  the  nearest-neighbor  distance  in 
AlgQj. 

For  theoretical  studies  the  (100)  surface  is  ideal:  calculated 
electronic  properties  can  be  compared  to  existing  Auger  and  UPS  results, 
and  calculated  bond  positions  and  lengths  can  be  compared  with  LEED, 
EXAFS,  and  EAPFS  results. 

Lang  and  Williams  (*0  have  studied  theoretically  the  adsorption  of 
oxygen  onto  a  jellium  surface.  While  jellium  is  a  good  approximation  to 
bulk  aluminum,  it  is  not  adequate  to  describe  the  local  adsorption 
process.  For  this  reason,  cluster  studies,  which  are  particularly  well 
suited  to  describe  such  local  phenomena  (  ),  have  been  performed. 
Harris  and  Painter  (  )  and  Messmer  and  Salahub  (  ),  have  studied  the 
adsorption  of  atomic  oxygen  onto  the  hole  site  of  the  A1  (100)  surface, 
reporting  electronic  structures  which  agree  with  existing  UPS  data  for 
the  initial  adsorption  of  oxygen  onto  this  alimunum  surface.  Studies  of 
the  adsorption  of  oxygen  at  the  remaining  two  high-symmetry  sites  of 
this  surface  have  not  previously  been  reported.  An  accepted  explanation 
of  the  incorporation  of  oxygen  atoms  beneath  the  aluminum  surface  dees 
not  yet  exist:  the  calculated  electronic  structure  and  binding 
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potential  curves  for  the  on-top  and  bridge  sites  presented  here  should 
shed  light  on  this  problem. 

In  this  study,  the  interaction  between  adsorbate  atoms  (in  this 
case  oxygen)  and  the  aluminum  (100)  surface  is  theoretically 
investigated  by  modeling  the  surface  with  a  small  number  of  atoms  and 
then  usir'3  the  unrestricted  Hartree-Fock  approximation  to  calculate  the 
electronic  structure  and  potential  energy  curves.  From  these  potential 
energy  curves  the  binding  energy,  bond  distance,  and  force  constants  of 
the  interaction  are  determined.  The  calculated  one-electron  orbitals 
give  added  insight  into  the  chemisorption  process,  and  calculated 
densities  of  states  are  compared  with  photoemission  data.  This 
calculational  procedure  is  able  to  accurately  describe  the  localized 
bonding  of  an  adsorbate  onto  a  surface  site,  the  aspect  of  chemisorption 
considered  most  important  (  ). 

The  localized  nature  of  the  chemisorptive  bond  is  indicated  by 
several  kinds  of  experimental  evidence.  An  adsorbed  atom  or  molecule 
can  bo  observed  to  hop  from  one  localized  site  to  another  (  ). 
Infrared  spectroscopy  studies  have  found  that  the  vibrational  spectra  of 
intermediates  on  surfaces  are  often  very  similar  to  the  spectra  of 
isolated  molecules  (  ).  These  experimental  results  reinforce  the 
intuitive  notion  that  the  chemisorptive  bond  is  similar  in  nature  to  the 
familiar,  localized  chemical  bond. 
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In  general,  the  loss  of  symmetry  in  the  direction  perpendicular  to 
the  surface  leads  to  the  choice  between  two  drastic  approximations  in 
theoretically  modeling  the  surface-adsorbate  interaction.  The 
maintenance  of  translational  symmetry  in  the  plane  parallel  to  the 
surface  greatly  simplifies  the  computational  difficulties  and  allows  for 
the  use  of  the  powerful  methods  developed  for  bulk  band  structure 
calculations.  The  alternative  approximation,  simulating  a  small  section 
of  the  surface  by  a  cluster  of  atoms,  can  more  effectively  describe  the 
localized  interaction  of  the  surface  with  a  single  atom. 

Numerous  calculations  of  the  electronic  structure  of  surfaces  using 
the  approximation  of  a  semi-infinite  surface  have  been  reported  (  ). 
The  major  drawback  of  this  approach  is  that,  in  general,  a  semi-infinite 
monolayer  of  adsorbed  atoms  must  be  considered  in  a  chemisorption 
calculation.  This  limitation  is  a  direct  result  of  the  symmetry 
assumed,  and  therefore  only  qualitative  calculations  can  treat  a  single 
adatom  on  a  semi-infinite  surface. 

By  relaxing  the  symmetry  requirements  of  the  semi-infinite  surface, 
one  gains  increased  flexibility  and  freedom.  For  instance,  atoms  can  be 
brought  down  over  a  variety  of  possible  bonding  sites.  Roughened  and 
stepped  surfaces,  exposed  corners,  and  small  particles,  which 
experimental  work  indicates  are  often  sites  favored  for  chemisorption, 
can  easily  be  simulated  by  the  cluster  approach.  Use  of  the 
semi-infinite  surface,  by  contrast,  is  limited  to  studies  of  the  perfect 
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Once  the  appropriate  cluster  has  been  chosen  to  siculate  the 
surface,  the  electronic  structure  of  the  cluster-adsorbate  systea  must 
be  calculated.  The  method  used  here  is  the  unrestricted  Hartree-Fock 
approximation,  froa  vrtiich  electronic  wavefuncbions,  electron  energy 
eigenfunctions,  and  the  total  electronic  energy  arc  calculated.  This 
calculation  is  repeated  at  varying  adsorbate-surface  separations,  and 
the  potential  energy  curve  can  be  generated  from  the  differences  in 
total  electronic  energy. 

The  Hartree-Fock  approximation  as  the  method  of  calculation  is 
attractive  due  to  two  important  advantages.  The  method  is  a  convenient 
first  step  towards  an  ab  initio  solution  of  the  many-particle 
Schrodinger  equation  for  a  fernion  system,  and  gives  a  useful 
zeroth-order  wavefunction  for  a  perturbation  calculation  of  the 
many-Lody  problem.  Second,  it  is  the  most  advanced  calculational  method 
in  which  elementary  physical  intuition  is  applicable.  In  more  complete 
calculations,  there  is  no  longer  a  one-to-one  correspondence  between 
particles  and  one-particle  wavef unctions,  and  the  independent  particle 
approximation  is  no  longer  applicable. 


1.2  TilS  ELECTRONIC  PROPERTIES  OF  THE  SILICA  SURFACE 


The  electronic  structure  of  crystalline  silica  (Si02)  ic  of 
technological  interest  in,  for  example,  the  manufacture  of  solid  state 
electronic  devices.  Consequently,  there  has  been  considerable 
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experimental  effort  in  this  area.  In  a  recent  review  of  experimental 
and  theoretical  results,  Griscoa  (  )  has  compiled  much  of  the  existing 
data  on  the  bull:  properties  of  silica.  The  band  gap  has  been  found  to 
be  about  8.9  eV,  with  a  valence  band  width  of  approximately  11  eV. 
Ultra-violet  photo  emission  spectroscopy  (UPS)  studies  by  Ibach  and  Rowe 
(  )  have  detected  no  occupied  surface  states  within  the  band  gap. 
Killians  (  )  has  studied  silica  surfaces  grown  on  crystalline  Si  and 
found  traps  at  2  eV  below  the  bulk  silica  conduction  band  edge.  Thes* 
traps  act  like  Coulomb  centers  with  a  positive  charge.  It  is  with  the 
siliv..  surface  and  its  defects  that  the  current  investigation  is 
concerned. 

Several  recent  experimental  studies  have  investigated  surface 
states  of  silica  (  ).  Schwidtal  (  )  has  found  a  radiation-induced 
feature  at  91  eV  on  the  high-energy  side  of  the  Si  L23  VV  (Vsvalence) 
transition  in  the  Auger  electron  spectrum  (AES),  and  suggests  that  it  is 
due  to  an  Si  L23  VD  (D=dcfect)  transition.  He  further  suggests  that 

f  *3 

this  defect  is  the  E3  center,  which  is  a  dar;gling  singly-occupied  spJ 
orbital  of  silicon.  Ibach  and  Rowe  (  ),  Fujiwara  and  Ogata  (  ), 
Lieske  and  Hezcl  (  ),  and  Bermudez  and  Ritz  (  )  have  measured  electron 
energy  loss  spectra  (ELS)  and  found  peaks  in  the  second  derivative 
spectrum  at  about  3.5,  5.0,  and  7.4  eV  (Fig.  1).  Ibach  and  Rowe  have 
attributed  these  transitions  to  a  partially  oxidized  surface  region, 
SiOx,  where  1<x<2.  Fujiwara  and  Ogata  (  )  have  shown  the  states 
associated  with  these  transitions  to  be  at  the  surface,  and  concluded 
they  were  due  to  metastable  SiOx.  Lieske  and  Hezel  (  )  have  associated 
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the  final  states  with  Si-Si  bonds  found  in  SiOx.  Host  recently, 
Bermudez  and  iiitz  (  )  have  studied  this  transition,  and  found  the  final 
state  to  be  due  to  "a  chemically  stable  species  formed  after  rupture  of 
the  silica  network,"  and  attributed  these  properties  to  a  surface  SisO 
double  bond.  For  this  bond  to  exist,  the  silica  surface  must  undergo 
reconstruction.  Evidence  for  such  reconstruction  from  low  energy 
electron  diffraction  (LEED)  measurements  is  due  to  Janossy  and  Menyhard 
(  ).  Also,  Hochstrasser  and  Antonini  (  )  have  observed  what  they 
believe  to  be  recombination  luminescence  due  to  the  rearranging  of  the 
surface  Si-0  bonds  immediately  after  cleaving.  They  measured  the 
lifetime  of  this  luminescence  to  be  less  than  10"®  seconds. 

Theoretical  studies  of  the  silica  bulk  electronic  properties  have 
been  performed  (  ),  and  the  calculated  densities  of  states  are  in 
good  agreement  with  experiment,  except  for  the  AES  peak  at  91  eV  and  the 
ELS  peaks  at  3.3,  5.0,  and  6.8  eV,  Both  of  these  discrepancies  can  be 
accounted  for  by  allowing  for  the  existence  of  islands  of  silicon  in  the 
bulk  Si02.  Bennett  and  Roth  (  )  have  calculated,  using  the  Huckel 
approximation,  the  electron  *  properties  of  many  clusters  which 
approximate  silica,  and  in  particular  gave  attention  to  defects.  These 
calculations  were  performed  in  the  bulk;  however,  the  clusters  were 
sufficiently  small  that  they  can  be  interpreted  in  terms  of  surface 
defects,  or  in  terms  of  a  surface  SiOx  region.  For  oxygen  deficient 
clusters  with  a  relaxation  of  the  silicon  atoms  near  the  defect,  the 
calculated  energy  levels  showed  transitions  at  approximately  the  same 
energy  values.  These  are,  however,  probably  due  to  small  clusters  (2-3 
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atoms)  of  silicon  in  the  SiOs.  Laughlin,  et  al .  (  ),  using  a  tight 
binding  method,  have  calculated  similar  defects  and  have  also  found 
states  which  appear  in  the  gap,  at  about  5  eV  above  the  valence  band 
edge.  They  have  also  studied  the  e'  center,  and  found  a  level  lying 
just  below  the  bulk  conduction  band.  Pollmann  an .  Pantelides  (  )  have 
performed  a  calculation  using  a  Green's  function  formulation  on  the 
ideal  terminated  cubic  (^-cristobalite)  surface.  Using  an  admittedly 
crude  racdel  for  the  surface,  they  found  no  states  in  the  band  gap. 

The  purpose  of  the  current  investigation  is  to  study  the  silica 
surface,  giving  particular  attention  to  the  three  ELS  peaks  observed  at 
3*5,  5,  and  7  cV.  The  cluster  model  is  again  employed,  for  the  same 
reasons  given  in  section  1.1.  Of  the  three  models  proposed  by 
experimentalists  to  account  for  these  transitions,  the  only  one  which 
has  previously  been  calculated  theoretically  is  theSiOx  model.  In 
chapter  4  are  presented  calculations  of  the  one-electron  energy  levels 
for  the  remaining  two  models,  the  center,  and  the  Si=0  surface 
double  bond.  A  comparison  is  then  made  between  the  three  models  and 


experiment. 


Bermudez  &  Ritz 


Energy  (eV) 


METHOD  OF  CALCULATION 


An  exact  solution  to  the  wave  equation  for  a  large  cluster  of  atoms 
is  not  at  present  attainable.  In  fact,  such  a  solution,  while  nice  in  a 
mathematical  sense,  would  be  difficult  to  interpret  physically.  The  ab 
initio  Hartree-Fock  theory  provides  an  approximate  method  of  solution  to 
determine  the  electronic  structure  of  a  large  cluster  of  atoms  which  is 
both  easily  handled  by  contemporary  computers  and  well  suited  for 
straightforward  interpretation. 


2.1  HARTREE-FOCK  THEORY 


Ab  initio  Hartree-Fock  theory  is  used  to  calculate  the  electronic 
wavefunctions  and  energies  for  the  finite  clusters  studied  in  this  work. 
The  solids  studied  contain  atoms  of  low  atomic  number;  therefore 
relativistic  effects  are  small  and  will  be  ignored.  Using  the 
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nonrelativistic  Schrodinger  theory,  the  many-body  wavefunction  x,  which 
is  a  function  of  electronic  and  nuclear  space-spin  coordinates,  is 
determined  by  the  equation 


where  H  is  the  raany-body  Hamiltonian 

H  :  T  ^,,,4 

with 


i! 


i. 


4* 


Here  we  have  used  atomic  units,  wherq^h',  me,  and  e=1 ,  the  unit  of  length 
is  the  bohr  (1  bohr 0.529A) ,  and  the  energy  is  given  in  Hartree  (1  Hy 
3S27.2  eV).  Upper  case  letters  refer  to  nuclear  properties:  Rj  is  the 
position  of  the  nucleus,  and  Mj  is  its  mass.  Lower  case  letters 
refer  to  electronic  properties:  is  the  coordinate  of  the  i^ 
electron.  The  four-vector  x^  denotes  both  the  spatial  coordinates  and 
the  spin  of  the  i^1  electron. 

The  Schrodinger  equation  is  simplified  by  employing  the 
Born-Oppenheiraer  approximation  (  )  to  separate  the  nuclear  and 
electronic  coordinates.  The  total  wavefunction  is  assumed  to  be 
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separable: 


t  a,h  - 


The  nuclear  wavefunction  'j'(X) ,  is  now  a  solution  of  the  equation 


j~~ T~%  ^ |  $&)  =  f|iJ) 


p-V 


while  the  electronic  wavefunction, '-t'(x,X),  is  a  solution  of  the  equation 
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A  solution  of  this  electronic  Schrodinger  equation  depends  on  the 
nuclear  coordinates  R  and  in  turn  defines  a  potential  energy  funotion  of 


the  coordinates: 


r.r 

which  deteraines  the  notion  of  the  nuclei.  The  problem  is  now  reduced 
to  solving  the  electronic  Schrodinger  equation  for  a  particular  set  of 
nuclear  coordinates.  This  equation  cannot  in  general  be  solved  exactly; 
one  therefore  takes  refuge  in  the  independent  particle  model  for  a 
method  of  obtaining  an  approximate  solution. 

This  approach  is  mathematically  equivalent  to  a  separation  of 
variables  technique.  The  raany-eloctron  wavefunction,  which  is  a 
functional  of  the  space-spin  coordinates  of  all  of  the  electrons,  is 
assumed  to  be  a  product  of  orbitals  which  are  functions  of  one-electron 
space-spin  coordinates.  However,  the  Pauli  principle  requires  that  the 
many-electron  wavefunction  be  antisymmetric  under  the  exchange  of  any 
two  electrons.  Therefore,  the  total  electronic  wavefunction  is 
approximated  by  an  antisymmetrized  product  of  one-electron  orbitals: 


where  $.is  the  antisyametrizer.  Since  the  Schrodinger  equation  is  atill 
not  in  general  solvable  analytically,  the  variational  principle  is 
introduced  to  obtain  approximate  solutions.  The  variational  principle 
states  that  the  normalized  expectation  value  of  the  Hamiltonian  using 
approximate  wavef unctions  is  a  rigorous  upper  bound  to  the  exact  energy. 
The  approximate  wavefunction  can  then  be  varied  to  minimize  the  energy 
and  thereby  be  the  best  approximation  to  the  exact  wavefunction. 
Applying  the  variational  principle  to  the  electronic  Schrodinger 
equation,  and  using  an  antisymaetrized  product  of  one-electron  orbitals, 
we  obtain  the  Hartree-Fock  equations: 


f  4  *  £ -l,-,  f  jw'J 
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vhe~e  F  is  the  Fock  operator: 
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The  matrix  elements  are  Lagrange  multipliers,  and  have  been 
introduced  to  insure  the  orthonormality  of  the  one-electron  orbitals. 
In  general,  a  unitary  transformation  can  be  performed  on  the  orbital 
space  which  diagonalizes  this  matrix.  The  Hartree-Fock  equations  can 
then  be  written  in  the  standard  form: 


I*-  »  •»»»*, 

1  )  ' 


These  integro-differential  equations  are  normally  solved  via  ar< 
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iterative  procedure.  Thus,  the  electron  orbitals  are  determined  in  the 
average  field  of  all  the  other  electrons.  The  only  constraints  imposed 
on  the  wavefunction  are  that  the  orbitals  be  orthonormal  and  be 
functions  of  the  space-spin  coordinates  of  only  one  electron.  This 
level  of  approximation  is  called  generalized  Hartree-Fock  (GKF)  and  has 
yet  to  be  solved,  so  further  constraints  are  imposed  to  simplify  the 
calculation.  The  simplest  constraint  which  can  be  imposed  is  to  require 
that  the  orbitals  be  eigenfunctions  of  Sz,  i.e.,  the  orbitals  be  of  the 
form: 


c  * 

k  (T) =  /  - 
1  ( f 


/2.I2,) 


where  ^(r)  is  the  spatial  part  of  the  orbital,  and  U  and  ^  are  the 
spin-up  and  spin-down  eigenfunctions  of  Sz,  respectively.  The  form  of 
the  Hartree-Fock  equations  is  unaffected  by  this  constraint,  and  this 
level  of  approximation  is  known  as  the  unrestricted  Hartree-Fock  (UHF). 
The  UHF  method  has  been  used  for  numerous  practical  calculations  of 
atoms,  molecules,  and  3olids  (  ),  and  is  the  method  used  in  this 


study. 


Additional  restrictions  may  be  placed  upon  the  one-electron 
orbitals  in  order  to  further  simplify  the  calculations.  The  spatial 
parts  of  the  orbitals  may  be  required  to  be  symmetry  eigenfunctions  of 
the  cluster  under  consideration,  i.e.,  each  orbital  must  transform  as  an 
irreducible  representation  of  the  space  group  of  the  cluster.  Finally, 
for  doubly  occupied  orbitals  the  spin-up  and  spin-down  spatial  orbitals 
may  be  assumed  to  be  equivalent.  This  level  of  approximation  i3  called 
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the  restricted  Hartree-Fock  method  (RHF).  One  advantage  of  RHF  is  that 
the  one-electron  orbitals  are  eigenfunctions  of  both  spatial  symmetry 
and  spin.  However,  for  singly  occupied  orbitals  complications  arise  in 
that  the  Fock  operator  differs  for  open  and  closed  shells  of  the  same 
symmetry,  thus  not  allowing  the  single  determinant  to  properly  describe 
the  wavefunction.  The  UHF  method,  while  not  necessarily  obtaining  a 
wavefunctinn  which  is  an  eigenfunction  of  spin  and  spatial  symmetry, 
does  give  a  lower  (i.e.,  better)  eigenfunction  of  energy.  Since  the 
"best"  solution  to  equation  2.6  is  the  wavefunction  which  gives  the 
lowest  energy,  the  UHF  method  has  therefore  been  used  in  this 
calculation. 


2.2  CLUSTER  MODEL 


Localized  phenomena  in  solids,  such  as  point  defects  or  surface 
adsorbates,  can  be  modeled  effectively  by  using  a  finite  cluster  of 
atoms  (  ).  Since  the  cluster  is  intended  to  represent  a  much  larger 
system,  one  must  apply  appropriate  boundary  conditions  for  the  cluster 
to  simulate  the  effects  of  the  environment.  These  have  been  discussed 
by  Kunz  and  Klein  (  ),  and  will  be  briefly  reviewed  here. 

The  solution  to  the  Hartree-Fock  equation  (2.11)  is  sought  for  a 
finite  cluster.  Suppose  the  region  of  the  cluster  is  termed  A,  and  the 
remainder  of  the  system,  the  environment  of  A,  is  termed  E.  The  problem 
is  to  partition  the  system  rigorously  into  a  cluster  and  an  environment. 
Let  us  consider  the  method  of  local  orbitals  of  Adams,  Gilbert,  and 
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Kunz(  ).  Here  the  Fock  operator,  F,  of  the  entire  system  is  divided 
into  a  part  which  includes  the  kinetic  energy,  nuclear  attraction  of  the 
electrons  and  nuclei  inside  A,  along  with  the  electron-electron 
potential  for  the  electrons  assigned  to  A,  F^;  and  another  part  which 
is  the  potential  in  A  due  to  the  environment,  U^: 


It  is  desired  to  study  only  part  of  the  system,  so  Instead  of  the  normal 
Hartree-Fock  equation  (2.11),  consider  instead  the  local  orbitals 
equation: 


where  W  is  an  arbitrary  Hernitian  operator.  Consider  now  an  ionic 
system.  The  potential  due  to  the  environment,  U^,  may  be  divided  into 

U 

two  parts:  V£  is  an  ionic  (Madelung)  contribution  and  is  long  range, 
and  is  the  remainder  and  i3  short  range.  Let  W  =  V^.  Making  use 
of  the  projector  properties  of  ,  and  solving  for  the  occupied 

orbitals,  equation  (2.14)  becomes 


Considering  only  the  orbitals  of  equation  (2.15)  which  lie  in  A,  the 
solutions  should  penetrate  only  weakly  into  E.  Because  does  not 
appreciably  penetrate  A,  and  because  ^  is  cancelled  by  ±  in  the 
limit  of  self-consistency,  the  approximate  equation  for  the  cluster, 
including  the  interaction  with  the  remainder  of  the  system,  is  then 
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This  approximate  equation  is  simply  the  equation  for  the  vavefunction  in 
the  cluster  of  interest:  imbedded  in  the  Madelur.g  field  of  the  remainder 
of  the  solid.  The  eigenvalues  "nj[  of  equation  (2.1*0  or  the  approximate 
equation  (2.16)  represent  the  eigenvalues  Cl  of  the  infinite  solid  only 
in  the  limit  that  the  orbitals  >j^  are  localized  in  the  region  A. 


A  basis  function  expansion  is  used  to  determine  the  unrestricted 
Hartree-Fock  orbitals.  Here  gaussian  type  orbitals  have  been  used  for 
this  expansion.  This  type  of  function  offers  the  advantage  of  having 
analytic  solutions  for  the  required  integrals;  however,  ujiny  such 
functions  are  required  to  achieve  the  desired  accuracy.  Cartesian 
gaussian  functions  have  been  used.  The  spatial  part  of  these  functions 
has  the  form: 
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The  determination  of  the  gaussian  exponents,  .rfj,  for  atoms  and  ions  has 
been  the  subject  of  considerable  study,  and  tabulated  sets  exist  in  the 
literature.  An  all-electron  calculation  including  the  core  electrons  of 
aluminum  and  silicon  would  be  impractical,  so  the  ab  initio  effective 
potential  of  Topiol,  et  al.  (  )  is  used  to  replace  these  core 
electrons.  The  integrals  are  performed  using  the  standard  POLYATOM 
integrals  progre  and  the  unrestricted  Hartree-Fock  calculations  are 
done  with  the  G.  T.  Surratt  program  UHFONE. 


THE  INITIAL  ADSORPTION  OF  OXYGEN  ONTO  THE 
ALUMINUM  (100)  SURFACE 


Aluminum  crystalizes  in  the  FCC  form,  with  a  lattice  constant  of 
*1.05  Aluminum  is  a  metal,  and  the  bulk  electronic  strunture  is 
approximated  very  well  by  the  wjelliumM  model.  However,  attempts  to 
employ  this  model  to  describe  the  local  chemisorption  process  (10)  have 
produced  results  which  are  not  consistent  with  experiment  (6).  Cluster 
techniques  have  been  used  to  study  the  chemisorption  process  (11, 12, 140) , 
with  results  which  agree  well  with  experiment  if  care  is  taken  in 
choosing  the  cluster. 

Chemisorption  is  a  localized  process,  and  as  such  the  cluster  model 
should  provide  an  accurate  description.  The  difficulties  that  one 
encounters  here  are  resultant  from  the  fact  that  aluminum  is  a  metal, 
and  as  such  is  not  well  described  by  a  small  cluster  of  atoms  (fig. 
3.1).  One  must,  therefore,  take  into  account  the  size  of  the  cluster 
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used  to  describe  the  substrate.  These  effects  will  be  discussed,  along 
with  the  binding  energy  of  the  oxygen  atom  for  each  surface  site 
considered.  Calculated  one-electron  energy  levels  for  tnese  binding 
sites  are  compared  with  the  results  of  photoemission  experiments. 


3.1  RESULTS  OF  CALCULATIONS 


Tne  positions  of  the  aluminum  centers  in  the  clusters  are  chosen  to 
reproduce  the  unreconstructed  structure  of  the  aluninun  (100)  surface. 
It  has  been  assumed  here  that  the  lattice  does  not  undergo  signifigant 
distortions  at  the  surface.  This  has  been  shown  to  be  true  from 
comparison  of  experimental  and  theoretical  LEED  studies  of  clean 
aluminum  surfaces  (M). 

As  possible  surface  sites  for  the  adsorbed  oxygen  atom,  considered 
here  are  the  three  high-symmetry  points  of  the  (100)  surface  (fig. 
3.2).  These  are  a  position  directly  above  an  aluminum  center  of  the 
first  substrate  layer  (on-top  position);  the  position  central  to  four 
aluminum  atoms  of  the  first  layer,  above  a  second  layer  atom  (hole 
position);  and  above  (or  below)  the  midpoint  of  a  line  connecting  two 
nearest  neighbor  aluminum  centers  (bridge  position). 

The  interaction  of  oxygen  with  the  aluminum  surface  at  the  on-top 
position  has  been  modeled  here  using  two  clusters:  A10,  with  a  single 
aluminum  atom  representing  the  surface;  and  AlgO,  with  five  first  layer 
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a to as  and  four  second  layer  atoms  modeling  the  surface.  The  binding 
energy  of  the  oxygen  atoa  to  the  aluminum  surface  is  shown  as  a  function 
of  its  distance  above  the  surface  in  figure  3*3*  For  both  surface 
clusters,  the  equilibrium  distance  is  about  3*5  bohr.  The  binding 
energies  differ  between  the  two  clusters;  this  energy  is  2.25  eV  for 
the  A10  cluster,  and  1.36  eV  for  the  AlgO  cluster.  The  charge  transfer 
from  the  aluminum  surface  to  the  oxygen  atom  is  about  0.7  electron  in 
both  cases. 

For  the  hole  site,  the  aluminum  surface  is  modeled  with  two 
clusters:  Al^o,  which  has  four  aluminum  centers  in  the  first  surface 
layer,  and  the  central  aluminum  in  the  second  layer;  and  AlgO,  which  is 
the  AI5O  cluster,  to  which  four  additional  second  layer  aluminum  centers 
are  added.  Tnc  importance  of  including  the  second  layer  atoms  at  this 
site  has  been  shown  by  previous  theoretical  studies  of  the  interaction 
of  oxygen  with  metal  surfaces  (45).  The  oxygen  binding  energies  as  a 
function  of  the  distance  above  the  surface  plane  are  shown  for  these 
clusters  in  figure  3*4.  Both  clusters  have  a  minimum  total  energy 
(maximum  binding  energy)  when  the  oxygen  center  is  about  0.2  bohr  below 
the  plane  of  the  first  surface  layer.  The  cluster  with  only  one  second 
layer  aluminum  atom  is  found  to  bind  the  oxygen  adsorbate  by  5 *4  eV, 
with  a  charge  transfer  to  the  oxygen  of  1.4  electrons.  The  larger 
cluster  is  found  to  bind  the  oxygen  by  4.8  eV,  with  a  charge  transfer  of 
1.3  electrons  to  the  adsorbed  oxygen. 
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“Hie  binding  energies  computed  for  this  site,  as  well  as  those  for 
the  bridge  situ,  must  be  carefully  extracted  from  the  calculated  total 
energies.  Since  the  oxygen  stow  penetrates  into  the  cluster  of  aluminum 
centers,  the  basis  functions  used  to  describe  the  electronic 
wave function?  associated  with  the  oxygen  adorbatc  may  inprove  the 
description  of  the  substrate  electronic  structure,  thus  lowering  its 
total  energy.  Tbi3  improvement  in  the  substrate  wavefuncticn  can  be 
integrated  into  the  caculation  of  the  binding  energy  by  simply  including 
the  basis  functions  of  tho  oxygen  atom  when  calculating  the  total  energy 
of  the  substrate  cluster. 

The  bridge  site  has  been  modeled  here  with  three  different 
clusters.  The  first  cluster  consists  of  two  nearest  neighbor  aluminum 
centers  of  the  surface,  and  an  oxygen  atom  above  the  midpoint  of  the 
line  connecting  these  centers.  The  second  cluster  contains  these  two 
centers,  as  well  as  two  additional  centers  beneath  them  from  the  second 
aluminum  layer.  For  the  final  cluster,  two  aluminum  centers  from  the 
third  layer  are  added.  The  inclusion  of  the  deeper  layer  aluminum 
centers  at  this  bonding  site  is  necessary  because  the  oxygen  atom 
penetrates  the  surface,  and  can  move  vertically  in  a  "tunnel"  between 
pairs  of  aluminum  centers  (see  fig.  3*2).  The  binding  energies  of  an 
oxygen  atom  interacting  with  these  aluminum  clusters  are  shown  in  figure 
3.5.  For  the  cluster  with  two  aluminum  centers,  the  oxygen  adsorbate  is 
bound  1.3  bohr  above  the  plane  of  the  surface.  A  charge  of  1.25 
electrons  is  transfered  to  the  oxygen  atom,  which  is  bound  by  3*1  eV. 
With  the  inclusion  of  the  two  second  layer  aluminum  centers  in  the 
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cluster,  the  equilibrium  position  of  the  oxygen  center  moves  to  "i  .65 
bohr  beneath  the  surface,  which  is  just  above  halfway  between  the  two 
aluminum  layers.  The  binding  energy  is  found  to  be  3*8  eV,  with  a 
charge  transfer  of  1.4  electrons  to  the  oxygen.  Inclusion  of  the  third 
aluminum  layer  produces  no  signifigant  change  in  these  results; 
however,  this  cluster  does  show  that  the  oxygen  center  does  not  want  to 
move  deeper  into  the  bulk  material,  as  the  total  energy  of  the  system  is 
found  to  be  lowest  for  an  oxvgen  center  position  just  below  the  first 
surface  layer. 


3.2  DISCUSSION  OF  RESULTS 


For  the  initial  stages  of  adsorption,  the  oxygen  adsorbate  is 
clearly  bound  below  the  aluninum  surface.  The  binding  is  strongest  at 
the  hole  site.  An  electric  charge  of  about  1.35  electrons  is  also 
transfered  to  the  adsorbed  oxygen.  These  results  are  in  agreement  with 
previously  reported  work  function  studies,  in  which  a  reported  decrease 
in  the  work  function  of  the  aluminum  surface  with  increasing  oxygen 
exposure  (up  to  one  monolayer  of  coverage)  has  been  attributed  to  the 
penetration  of  negative  oxygen  ions  beneath  the  aluminum  surface  (6). 

The  ultra-violet  photoemission  spectrum  has  been  reported  for  low 
oxygen  coverage  on  the  aluminum  (100)  surface  (5);  this  spectrum  is 
shown  in  figure  3*6.  Also  shown  is  the  projected  density  of  states  for 
the  oxygen  adsorbate  previously  calculated  by  Messmer  and  Salahub  (12) 
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for  an  oxygen  position  at  the  aluminum  surface  at  the  hole  site.  These 
authors  did  not  calculate  the  electronic  structure  for  an  oxygen 
position  below  the  surface,  nor  did  they  study  the  interaction  of  an 
oxygen  adsorbate  with  any  of  the  other  high-symnetry  (100)  surface 
sites. 

As  reported  elsewhere  (*10),  the  one-electron  energy  levels  of  this 
calculation  are  in  agreement  with  experiment  and  previous  theory.  These 
energy  levels  are  shown  in  figure  3.6  for  both  the  hole  and  bridge 
sites,  with  the  oxygen  position  corresponding  to  the  minimum  total 
energy  for  the  cluster.  For  both  sites,  the  energy  levels  reported  here 
are  those  of  the  largest  cluster  calculated.. 

Photoenission  experiments  studying  tho  aluminum  core  levels  at  low 
oxygen  exposures  have  been  reported,  with  a  reported  shift  in  the 
aluminum  2p  energy  of  1.3eV  toward  lower  energy  (5).  In  the  present 
calculation,  the  shift  of  the  alurainn  core  levels  is  found  to  be  1.1  eV 
in  the  same  direction.  Tne  cause  of  this  core  level  shift  is  the 
electric  dipole  produced  by  the  charge  transfer  from  the  aluminum 
substrate  to  the  oxygen  atom  upon  adsorption.  It  has  been  noted  by 
Flodstrom,  et  al.  that  the  aluminum  core  level  shift  increases  to  a 
value  of  2.6  eV  at  higher  oxygen  exposure  (about  one  monolayer  of 
coverage),  thus  lowering  these  orbital  energies  to  the  values  found  in 
AlgOj.  These  results  are  evidence  for  a  two  stage  oxidation  process, 
with  the  first  step  being  chemisorption  at  the  hole  site,  and  the  second 
step  being  incorporation  of  oxygen  beneath  the  aluminum  surface. 
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The  agreement  between  the  calculated  one-electron  energies  for  the 
bole  site  presented  here,  along  with  those  of  Messmer  and  Salahub,  and 
reported  experimental  UPS  spectra  tend  to  support  the  idea  that  the 
initial  chemisorption  of  oxygen  occurs  at  the  hole  site,  as  one  would 
expect  from  the  binding  energies  reported  in  section  3*1* 

The  hole  site  is  also  favored  by  extended  appearance  potential  fine 
structure  (BAPFS)  measurements.  Using  this  technique,  den  Boer,  et  al. 
(4)  have  reported  a  nearest  neighbor  oxygen-aluminum  distance  of  1.98  A, 
or  3.7*1  bohr.  Messmer  and  Salahub  have  reported  a  value  of  2.02  A  (3*82 
bohr),  and  the  results  of  the  calculations  presented  here  place  the 
oxygen  at  1.92  A  (3.63  bohr)  from  the  nearest  aluminum  center  for  the 
hole  site.  For  higher  oxygen  exposure,  den  Boer^et  al.  have  reported 
an  experimental  value  of  this  nearest  neighbor  distance  of  1.91  A*  This 
figure  is  probably  an  average  of  the  two  nearest  neighbor  spacings  found 
in  AlgOj;  these  being  1.86  A  (3*51  bohr)  and  1.97  A  (3*72  bohr).  The 
longer  of  these  spacings  corresponds  to  oxygen  binding  at  the  hole  site, 
while  the  shorter  corresponds  to  oxygen  adsorption  at  the  slightly  less 
energetically  favorable  bridge  site.  The  present  calculation  places  the 
nearest  neighbor  distance  at  the  bridge  site  at  1.71  X,  or  3*25  bohr. 

The  results  presented  above,  along  with  previous  theoretical  and 
experimental  results,  show  that  for  the  initial  interaction  of  oxygen 
with  the  aluminum  (100)  surface,  chemisorption  takes  place  at  the  hole 
site.  The  present  studies  have  extended  previous  theoretical 
understanding  of  the  interaction  of  oxygen  with  this  surface  to  include 
the  two  additional  high-symetry  (100)  sites.  Work  function,  UPS,  and 
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EAPFS  studies  have  reported  that  the  nature  of  this  interaction  changes 
when  the  oxygen  coverage  approaches  one  nonolayer;  the  theoretical 
results  presented  here  lead  one  to  conclude  that  at  this  coverage, 
oxygen  begins  to  adsorb  at  the  bridge  site,  since  the  cnergitically  more 
favorable  hole  sites  arc  already  occupied.  The  chemisorption  of  oxygen 
at  the  bridge  site  allows  for  the  incorporation  of  the  oxygen  adatoas 
beneath  the  aluminum  surface  at  high  levels  of  oxygen  exposure,  leading 
to  the  formation  of  the  oxide  AljOj. 


THE  ELECRONIC  PROPERTIES  OF  THE  SILICA  SURFACE 


Silica  occurs  in  six  different  crystalline  forms,  as  well  as  the 
amorphous  form.  All  forms  of  Si02  are  based  on  the  SiOjj  tetrahedron, 
and  differ  in  the  Si-O-Si  bonding  angle.  The  crystalline  form  most 
often  studied  experimentally  is  e(-quartz,  and  it  is  this  form  which  we 
have  used  here,  *£-quartz  is  hexagonal  in  structure,  with  three  Si02 
molecules  In  each  unit  cell.  The  Si-O-Si  bond  angle  is  IW* . 

The  silica  surface  may  be  described  accurately  using  a  cluster  of  a 
small  number  of  atoms.  Indeed,  calculations  on  the  bulk  material  (  ) 
have  used  the  cluster  approach.  To  take  into  account  the  effect  of  the 
neighboring  bulk,  appropriate  boundary  conditions  are  imposed,  as  has 
been  discussed  in  chapter  2.  In  the  case  of  silica,  the  cluster  is 
embedded  in  a  point  charge  array,  with  the  charge  transfer  determined  in 
a  self-consistent  manner.  This  gives  the  proper  electrostatic  potential 
in  the  region  of  the  surface  cluster,  and  also  provides  for  charge 
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neutrality.  Previous  calculations  have  cither  employed  a  free  space 
termination,  or  imposed  periodic  boundary  conditions  which  require 
defects  to  be  present  in  all  unit  cells  of  the  crystal. 

The  clusters  used  here  for  the  bulk  calculations  were  an  SiOq 
tetrahedron,  and  an  Si-O-Si  molecule.  Both  clusters  were  embedded  in  a 
point  charge  2rray  of  3*3*3  unit  cells  in  size.  The  calculations  of  the 
surface  were  done  using  a  cluster  consisting  of  the  surface  silicon 
atom,  its  two  neighboring  oxygen  atoms,  and  the  two  adjacent  silicon 
atoms  to  form  two  joined  Si-O-Si  molecules.  This  cluster  was  embedded 
in  a  point  charge  array  of  3*3*2  unit  cells.  The  position  of  the 
surface  silicon  atom  was  adjusted  in  order  to  minimize  the  total  energy. 
An  oxygen  atom  was  then  placed  above  the  surface  silicon  in  order  to 
study  the  surface  bonding  state. 

The  one-electron  energy  levels  calculated  for  the  bulk  SiC^  arc  in 
good  agreement  with  experiment  and  previous  theory  (fig.  4.1).  The 
calculated  band  gap  is  9.4  eV,  as  compared  to  8.9  eV  for  experiment. 
The  valence  band  width  is  calculated  to  be  9*2  eV,  while  experimental 
measurements  give  about  11  eV.  The  oxygen  2s  band  is  found  at  30.3  eV 
below  the  conduction  band,  while  experiment  places  it  28  eV  below.  The 
valence  bands  are  found  to  be  mostly  oxygen  2p,  with  some  silicon  3s  and 
3p  character  mixed  in.  This  is  expected,  and  agrees  with  the  charge 
transfer  of  nearly  one-half  of  an  electron  from  each  silicon  to  each  of 
4  oxygen  atoms,  or  a  configuration  of  Si2+o~,  as  determined  by  a 
Kulliken  population  analysis.  This  agrees  with  the  observation  that 
silica  has  both  ionic  and  covalent  properties. 
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For  the  unreconstructed  surface,  the  E'  center  is  found  to  have  a 

3 

singly  occupied  energy  level  1.7  eV  above  the  top  of  the  bulk  valence 
band  (fig.  4.2).  Stephenson  and  Binkowski  (  )  have  observed,  using 

XPS,  an  occupied  energy  level  at  0.75  eV  above  the  bulk  valence  band 
edge  for  samples  cleaved  in  vacuum.  They  believed  this  level  to  be  an 
intrinsic  feature  of  the  bulk  SiOg.  However,  their  sampling  depth  is 
only  about  30  A,  so  it  is  probable  that  this  level  is  an  E'  center  at  or 
near  the  surface  (  ).  The  method  of  sample  preparation  which  they  have 

used  (grinding)  could  easily  have  caused  these  broken  bonds  to  be  formed 
(  ).  It  is  energetically  favorable  for  an  oxygen  atom  to  bind  to  this 

surface  silicon  atom,  with  a  binding  energy  of  5.1  eV,  which  is 
considerably  greater  than  the  2.6  eV  per  oxygen  atom  necessary  to 
dissociate  molecular  oxygen.  We  shall  now  turn  our  attention  to  this 
configuration,  which  is  the  Si=0  surface  double  bond. 

The  occupied  one-electron  energy  levels  of  the  surface  state  are 
all  below  the  top  of  the  bulk  valence  band.  This  is  consistent  with 
optical  studies,  which  have  found  no  occupied  surface  states  in  the  band 
gap  (  ).  The  lowering  of  the  valence  band  can  be  attributed  to  the 

O-Si-O  bond  angle  at  the  surface  being  greater  than  the  perfect 

0 

tetrahedral  109  angle.  The  bottom  of  the  conduction  band  for  the 
surface  is  found  to  be  4.6  eV  below  the  bulk  conduction  band  edge.  The 
occupied  valence  surface  states  can  be  described  as  bonding  and 
non-bonding  states  between  the  adsorbed  oxygen  and  the  surface  silicon 
atom  (fig.  4.3a,b).  When  a  surface  electron  is  excited  out  of  the 
valence  band,  an  anti-bonding  orbital  (fig.  4.3c)  drops  out  of  the 
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conduction  band  and  is  singly  occupied.  Its  energy  lowers  considerably 
due  to  relaxation,  and  is  found  at  1.4  eV  above  the  top  of  the  bulk 
valence  band  (fig.  4.2).  The  position  of  the  oxygen  atom  is  0.1  a.  u. 
farther  out  from  the  surface  silicon  atom  for  this  excited  state, 
relative  to  its  ground  3tate  position.  This  energy  difference  is  in 
good  agreement  with  Bermudez  and  Ritz  (  ),  and  the  transition  from  a 
non-bonding  to  an  anti-bonding  orbital  of  the  surface  bond  is  as  they 
have  described.  The  energy  change  associated  with  this  transition  is 
6.1  eV,  which  is  in  agreement  with  the  experimental  value  of  5.0  eV  if 
correlation  effects  are  considered.  The  7*4  eV  transition  can  be 
understood  in  terms  of  valence  band  structure,  with  transitions  from 
levels  in  the  bonding  part  of  the  valence  band  to  the  same  anti-bonding 
orbital. 

It  is  noted  that  this  model,  the  Si=0  surface  bond,  does  not 
predict  the  surface  electronic  transition  at  3*5  eV  that  is  seen  in  the 
second  derivative  ELS  spectra.  It  has  been  shown  by  Gallon  and 
Underwood  (  )  that  this  peak  is  in  fact  an  artifact  of  the  second 
derivative  mode  of  detection,  and  it  i3  not  seen  in  the  non-derivative 
ELS  spectrum.  They  argue  that  the  3*5  eV  peak  is  produced  by  the 
overlap  of  the  "wings"  of  the  primary  and  the  5.0  eV  peaks.  The 
calculations  presented  here  agree  with  this  interpretation. 

One  remaining  point  of  controversy  is  the  difference  between  XPS 
and  UPS  in  discribing  the  top  of  the  valence  bands  of  Si02  (fig.  4. 4) 
There  exists  at  this  time  no  agreement  as  to  the  cause  of  this 
discrepancy.  The  calculations  presented  here  do,  however,  favor  one 
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explanation.  The  two  spectra  are  spatially  sensitive  to  different 
parts  of  the  sample:  XPS  is  sensitive  to  the  bulk  for  initial  states 
with  small  binding  energy,  while  IJ0.8  eV  UPS  is  surface  sensitive.  In 
both  the  XPS  and  UPS  techniques,  the  kinetic  energy  of  the  photoelectron 
is  sufficiently  large  that  the  final  state  may  be  approximated  by  a 
plane  wave.  Therefore  the  difference  between  the  XPS  and  the  UPS 
spectra  must  reflect  the  difference  between  the  bulk  and  surface  valence 
band  densities  of  states.  Comparison  between  the  XPS  spectrum  and  the 
calculated  bulk  density  of  states,  and  between  the  UPS  spectrum  and  the 
calculated  surface  density  of  states  (for  the  reconstructed  surface) 
shows  that  3uch  an  explanation  docs  indeed  account  for  the  difference  in 
the  two  spectra  (fig.  *4 . 4  > . 

The  reconstructed  silica  surface  has  been  studied  in  order  to 
determine  the  nature  of  the  low  energy  loss  ELS  peaks.  We  have  found 
agreement  between  the  theoretically  calculated  one-electron  energy 
levels  presented  here  and  previous  experimental  results.  Of  the  three 
models  which  can  account  for  the  peaks,  namely  a  partly  oxidized  surface 
region  SiOx,  the  center,  and  the  Si=0  double  bond,  we  conclude  that 
the  last  is  most  likely.  Support  for  this  conclusion  comes  from 
experimental  evidence  for  reconstruction,  and  from  the  large  binding 
energy  of  the  oxygen  atom  to  the  surface. 
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CONCLUSIONS 


In  this  study,  the  ab  initio  unrestricted  hartree-Fock  nethod  has 
been  used  to  calculate  the  bulk  and  surface  electronic  structure  of  of 
silica,  and  to  study  the  initial  interaction  of  oxygen  with  the  aluminum 
(100)  surface.  In  both  of  these  materials,  the  cluster  model  has  been 
shown  to  be  a  useful  tool  in  the  calculation  of  the  electronic 
structures. 

In  studying  the  interaction  af  oxygen  with  the  aluminum  (100) 
surface,  the  size  of  the  aluminum  cluster  used  had  a  considerable  effect 
on  the  results.  Since  aluminum  is  a  metal,  a  fairly  large  cluster  of 
atoms  is  needed  to  describe  the  substrate  accurately.  That  the  cluster 
model  worked  at  all  for  this  calculation  is  due  to  the  fact  that 
chemisorption  is  a  local  process. 
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It  has  been  shown  that  an  oxygen  adsorbate  is  bound  most  strongly 
at  the  hole  site  of  the  clean  aluninun  surface.  Tne  equilibrium 
position  of  the  adsorbed  atom  is  about  0.2  bohr  below  the  plane  of  the 
surface.  A  charge  of  1.3  electrons  is  transfered  from  the  aluminum 
substrate  to  the  oxygen  adsorbate,  Thi*  charge  transfer  creates  an 
electric  dipole  which  causes  the  aluminum  core  levels  to  be  shifted 
downward  in  energy  by  1.1  eV.  Previous  experiments  studying  the  work 
function,  UPS  spectrum,  and  the  EAPFS  have  reported  results  for  the 
initial  adsorption  of  oxygen  onto  this  surface  which  are  in  agreement 
with  the  results  of  the  calculations  presented  here. 

These  experimental  results  have  also  shown  that  the  nature  of  the 
adsorption  process  changes  at  about  one  monolayer  of  coverage.  The 
theoretical  studies  reported  here  lead  one  to  conclude  that  at  this 
coverage  the  energetically  most  favored  sites  for  chemisorption,  the 
hole  sites,  become  fully  occupied,  and  adsorption  continues  at  the 
bridge  site.  It  is  adsorption  of  oxygen  at  thi3  site  which  allows  for 
incorporation  of  oxygen  beneath  the  aluminum  surface,  and  leads  to  the 
formation  of  the  oxide  AlgOj. 

Also  reported  here  is  a  study  of  the  bulk  and  surface  electronic 
structure  of  silica  (/-quartz).  The  cluster  model  was  again  used.  The 
bulk  solid  and  the  surface  were  both  modeled  with  clusters  of  a  few 
atoms,  plus  the  appropriate  boundary  conditions.  In  this  case,  the 
boundary  condition  imposed  was  a  point  ion  array  to  provide  for  charge 
neutrality  and  to  provide  the  correct  Madelung  field. 
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For  the  silica  bulk,  the  calculated  one-elcctron  energy  levels  were 
compared  with  experimental  XPS  spectra,  and  the  agreement  was  found  to 
be  excellent,  as  reported  in  chapter  For  the  reconstructed  silica 
surface,  the  structure  of  the  Si=0  double  bond  was  found  to  be  the 
energetically  most  favorable  surface  configuration.  The  electronic 
structure  of  this  surface  bond  is  able  to  explain  the  low  energy  ELS 
peaks  seen  experimentally.  This  calculations  reports  a  value  of  6.1  eV 
for  the  lowest  energy  transition,  as  compared  to  5.0  eV  for  the 
experimental  result. 

Finally,  the  difference  in  the  reported  valence  band  structures  as 
measured  in  XPS  and  UPS  experiments  is  explained.  The  XPS  method  is 
sensitive  to  the  bulk  material,  while  UPS  is  surface  sensitive. 
Comparison  of  these  spectra  with  the  calculated  valence  band  energy 
levels  shows  good  agreement  between  the  theory  presented  here  and  the 
experimenal  spectra. 

Tne  ab  initio  unrestricted  Hartree-Fock  method,  along  with  the 
cluster  model,  has  been  shown  to  accurately  describe  the  local 
electronic  properties  of  many  systems.  In  this  report,  this  method  has 
been  shown  to  be  successful  in  describing  the  adsorption  of  oxygen  onto 
the  aluminum  surface,  and  in  describing  the  bulk  and  surface  electronic 
structure  of  SiC^.  This  technique  has  also  been  used  to  successfully 
calculate  the  electronic  structure  of  semiconducting  polymers,  including 
defects  and  impurities  ( h6 ) ,  and  in  the  description  of  the  excitonic 
structure  of  crystalline  silicon  (^7). 
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